Faraday effect in graphene enclosed in an optical cavity and the equation of motion 
method for the study of magneto-optical transport in solids 
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We show that by enclosing graphene in an optical cavity, giant Faraday rotations in the infrared 
regime are generated and measurable Faraday rotation angles in the visible range become possible. 
Explicit expressions for the Hall steps of the Faraday rotation angle are given for relevant regimes. 
In the context of this problem we develop an equation of motion (EOM) method for calculation of 
the magneto-optical properties of metals and semiconductors. It is shown that properly regularized 
EOM solutions are fully equivalent to the Kubo formula. 
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I. INTRODUCTION 

Electromagnetic radiation emitted by far stellar ob- 
jects travels for long periods of time through very di- 
luted concentrations of interstellar gases, traversing re- 
gions where week magnetic fields exist. In this circum- 
stance, the polarization of the electric field rotates due to 
its interaction with the gases immersed in the magnetic 
field. Due to the enormous traveling distances through 
such interstellar regions, the degree of rotation of the po- 
larization can be important. This magnetic rotational 
effect turns out to be a problem in astrophysics, since it 
modifies, in an unpredictable way, the polarization state 
of the emitted radiation, introducing additional difficul- 
ties in the interpretation of astronomical observations. 
In the electrodynamics of metals and insulators the ef- 
fect of polarization rotation induced by a magnetic field 
was first discussed by Faraday 1 and, on Earth, has many 
different applications. 

In magneto-optics, the effect coined optical Faraday 
rotatiorP refers to the rotation of the plane of polariza- 
tion of light when it transverses either a dielectricP or a 
metal,'^in the presence of a static magnetic field applied 
along the direction of propagation of the electromagnetic 
wave. In addition to the rotation of the plane of polar- 
ization, the polarization itself acquires a certain degree 
of ellipticity. In dielectrics, the effect can be explained 
using a model of harmonic oscillators coupled to light 
In metals, the effect has its roots in the Hall effect .31 

For a two-dimensional (2D) metal, such as graphene, 
in the Hall regime, the conductivity becomes a tensor 
a, with finite (nonzero) values for both diagonal and 
off-diagonal components. In magneto-optics, the com- 
ponents of the tensor depend both on the frequency of 
the impinging electromagnetic wave and on the cyclotron 
frequency of the electrons, due to the magnetic field per- 
pendicular to the plane of the metal. The response of the 
electrons to the external magnetic field has two regimes: 
(i) the semiclassical limit, of low fields and/or a high 



electronic density; and (ii) the quantum Hall regime, of 
strong fields and/or a low electronic density. 

For interpretation of the optical Faraday rotation, in 
the semi-classical regime, the Drude theory of metals 
suffices."^ In the case of graphene, it is possible to change 
its electronic density either by use of a gate or by the ad- 
sorption of molecules.'^ At high doping, graphene is in 
the semiclassical regime and Boltzmann transport theory 
can be used to compute the Hall conductivity.'^ 

In the absence of disorder and other relaxation mecha- 
nisms (such as electron-phonon scattering), the conduc- 
tivity of graphene (at zero magnetic field) would be exclu- 
sively determined by interband transitions. In the limit of 
no disorder, the optical conductivity of doped graphene, 
in the infrared region of the spectrum and at zero mag- 
netic field, is given bjJsHUl 

^ agHF^huj ~ 2Ef) , (1) 

where CTg = 7re^/(2/t) is t he so-called ac universal con- 
ductivity of graphene .ISliSHUl 

When a magnetic field is applied perpendicularly to 
graphene's surface, the system develops a finite Hall con- 
ductivity. In the quantum regime, it was shown that the 
Faraday rotation angle 6p is solely determined by the fine 
structure constant a, and presents a step-like structure as 
the Fermi energy crosses different Landau levels (LLs).l^ 
The estimated Faraday rotation steps' height in this case 
is of the order of 9f ~ 0.4° a magnitude that can be 
resolved experimentally.lSSl In the context of topological 
insulators, similar quantization rules in certain thin-film 
geometries have been derived in Refs. El] and EH We 
note in passing that, when the external magnetic field 
is absent, a dynamic Hall effect can still be induced by 
using circularly polarized light impinging on graphene at 
a finite angle with the normal to the graphene surface 

On the theoretical side, the magneto-optical trans- 
port properties of graphene ha ve b een investigated with 
the Green's function methocP*^, and by means of 
numerical implementations of the Kubo formula, us- 
ing exact diagonalization^^ and Chebyshev polynomial 
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expansions.'- - These approaches come with pros and cons: 
numerical studies allow the exploration of general scenar- 
ios, whereas Green's functions allows one to obtain ana- 
lytic results, but many times at the expense of lengthy 
calculations. 

Motivated by the need for analytical flexible analytical 
tools, the equation of motion (EOM) method employed 
in Ref . [35] is generalized to include the effect of a mag- 
netic field. As shown later, starting from a small set of 
EOMs, an adequate treatment permits the derivation of 
response functions with correct analytical properties (i.e., 
satisfying Kramers-Kronig causality relations). 

The present paper is divided into two main parts. In 
Sec.jnjwe present the EOM method for calculation of the 
magneto-optical transport in metals and semiconductors; 
to be concrete, the method is described in the context 
of the properties of graphene. In Sec. |III| we describe 
in detail the Faraday effect in graphene and propose an 
experimental setup that is able to enhance the Faraday 
effect up to the visible range. Section [Til] relies heavily 
on the results derived in Sec. |IT1 Some technical details 
are given in the Appendixes. 

We have chosen to organize the subjects according 
to the following interests of different readers: a reader 
having a primary interest in the Faraday effect, and fa- 
miliar with the details on the magneto-optical proper- 
ties of graphene, should be able to read Sec. |III| with a 
bird's-eye reading of Sec. |llj A reader interested in the 
Faraday effect in graphene but not well acquainted with 
its magneto-optical properties may want to go through 
Sec. first. Finally, reading Sec. |TT] alone may appeal to 
readers interested in applying the EOM method to an- 
other problem of interest bearing no relation to graphene. 



II. EQUATION OF MOTION METHOD FOR 
CALCULATION OF THE MAGNETO-OPTICAL 
CONDUCTIVITY 

Here, we develop the EOM approach to the calculation 
of the magneto-optical properties of a semiconductor. To 
be concrete, the method is presented in the context of the 
optical response of graphene. 

Electrons constrained to two dimensions are responsi- 
ble for a variety of quantum manifestations, a striking 
example being the integer quantum Hall effect (IQHE). 
Measured in semiconductor 2D electron gases more than 
30 years agcPSl and in the early days of graphene, in 
both monolayei'^1271 bilayer samples^ (very recently 
also in trilayer graphene^^), the static quantum Hall ef- 
fect is a hallmark of elementary excitations in electronic 
systems.!^ 

Its dynamical analog — the ac quantum Hall effect — 
can provide additional information about charge carriers, 
such as the opening of gaps in the spectrum .I^H Recent ad- 
vances in time-domain spectroscopy in the Thz regim^^ 
have paved the way to measurement of dynamical op- 
tical conductivities at impinging field energies closer to 



the scale of interest. The goal is to reach cyclotronic 
energies, usually 0(10) meV in fields of 1-10 T, where 
strong optical responses take place. The so-called optical 
quantum Hall conductivity of 2D electron gases shows a 
robust plateaux as the Fermi energy is swept, although no 
quantization rule for the plateaux's height exists.'^f Due 
to its peculiar band structure, graphene has been pre- 
dicted to display a characteristic optical quantum Hall 
effect which should be detectable via Faraday rotation 
measurements.^^ In the semiclassical regime, on the other 
hand, the Faraday rotation of graphene was reported to 
be 0(1) degrees in fields of a few teslaj^ a surprisingly 
high value for a one-atom-thick electronic system. 




FIG. 1: Lattice structure and Brillouin zone of monolayer 
graphene. Left: Hexagonal lattice of graphene, with the next 
nearest neighbor, Si, and the primitive, a^, vectors depicted. 
The area of the primitive cell is Ac = 3\/3ao/2 — 5.1 A^, 
and ao — 1.4 A. Right: Brillouin zone of graphene, with the 
Dirac points K and K' indicated. Close to these points, the 
dispersion of graphene is conical and the density of states is 
proportional to the absolute value of the energy. 



A. Graphene 

The starting point of the present analysis is the low- 
energy continuum description of single-layer graphene; 
having two (carbon) atoms per unit cell and sixfold sym- 
metry, its elementary excitations obey a 2D Dirac equa- 
tion with linear electronic dispersion.'^ This section is 
meant to fix the notation. The Brillouin zone of graphene 
has six corners, and among these, only two are inequiv- 
alent, the so-called K and K' Dirac points (see Fig. [l]). 
At these points, the valence and conduction bands touch, 
with a linear electronic spectrum up to energies of ^2 eV. 

We assume, in what follows, that the two Dirac points 
can be treated independently, and introduce the val- 
ley degeneracy index, = 2, when pertinent. This 
consideration is justified for typical experimental condi- 
tions (i.e., low concentrations of scattering centers, finite 
temperatures, etc.) and provides an accurate descrip- 
tion of gr aphe ne's electronic transport properties at finite 
densities.l ^^ l '^^ l 

In accordance, we resort to the 2x2 Dirac Hamilto- 
nian of graphene, describing the physics of elementary 
excitations within the K valley, Hj^ = vpcr ■ p, where 
vp — lO^m/s is the Fermi velocity, cr = (ax,cry) [with 
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ai {i = x,y) denoting Pauli matrices], and p is the mo- 
mentum of the low-energy excitation (measured relative 
to the K point) Hj^ has eigenvalues given by 



E = ±hvF\k\ , 



(2) 



[with k = {kx:,ky) denoting a 2D wave vector], and (nor- 
malized) wave functions given by 



V'A,k(r) 



1 

7n 



1 



(3) 



where A is the area of the graphene sample, A = 
-|-1(— 1) for electron(hole)-like excitation, and 6'k ~ 
arcta.n{ky / kx) ■ 

The electromagnetic field can be incorporated via min- 
imal coupling, p — p + eAg, where — e < is the 
electron charge, and the vector potential Ag relates to 
the electromagnetic field according to the usual relations, 
B = V X Ag and E = -dA^/dt. 

Here, the vector potential contains the information 
about the impinging electromagnetic radiation, and pos- 
sible external static magnetic fields. Assuming light lin- 
early polarized along the x axis, the radiation term reads 
A = [^o(r)e~"^* -I- c.c.Jex, where w stands for the fre- 
quency of the radiation field and ^o(r) describes its po- 
sition dependence. For clarity of exposition, we separate 
the light-matter interaction term from the free Hamilto- 



H = Hq + evpcr ■ A , 



(4) 



where Hq = + evpcr ■ As, with A^ describing the 
static magnetic field. 

A typical experimental scenario corresponds to a con- 
stant magnetic field B > applied in the transverse di- 
rection with respect to the graphene plane. In such case, 
LLs develop and the eigenenergies of charge carriers be- 
come quantized according tcpSl 



En = sign(n)- 



,n = 0,±1,±2,. 



(5) 



with Ib = (eB) denoting the magnetic length. 

Choosing the gauge Ab = {0,Bx,0) results in the fol- 
lowing set of Landau eigenfunctions. 



C„ 



(t>\n\-l{x) 

zsign(n)(/)|„|(a;) 



Jkyy 



(6) 



where (/.„(x) = e~^<-^^"/^H„{^{x))/^nl2'-^lB, Hn{x) is 
the Hermite polynomial of degree n > 0, (j)_i(x) = 0, and 
^(x) stands for the dimensionless center of the Landau 
orbit, ^{x) — iBky+x/ls- Here, L is the linear dimension 
of the system in the y direction and C„ is a normalization 
constant that distinguishes the zero-energy level from the 
remaining levels, Cn = 1 for n = and C„ = 1/a/2 for 
|n| > 1. 

Having reviewed the basics of the graphene's electronic 
low-energy theory, in what follows we present the EOM 
approach to the study of magneto-optical transport. 



B. Theoretical methods 



In the context of electronic systems, the EOM was 
extensively used in calculations of light polarization in 
semiconductor laser theory.!^ Recently, it has been used 
to study excitons in graphene in zero field. 

The EOM approach avoids the calculation of current- 
current correlators (i.e., Kubo formula), and, hence, pro- 
vides a shortcut to determination of the response of elec- 
tronic systems to external perturbations. As shown in 
detail in Appendix [C] with an appropriate regularization 
procedure, the EOM solutions become fully equivalent 
to the Kubo formula, and hence provide an accurate de- 
scription of transport in the linear response regime. An- 
other advantage of the present approach is that it allows 
for the calculation of non-linear corrections to the con- 
ductivity. 

At the heart of the EOM approach to calculation of 
the magneto-optical conductivity is the Heisenberg equa- 
tion for the electronic current density, J{t), in the pres- 
ence of an external electromagnetic field, i.e., dJ/dt — 
{i/h)[H, J], with H being the total Hamiltonian, Eq. Q. 
Having solved for the current density of the system in the 
presence of external perturbation, in first order in the ex- 
ternal field A, the optical conductivity follows from the 
constitutive electromagnetic relation 



(Jij{uj) = gs 



X — 



(7) 



where 0(w) relates to the average 0{t) [O — Ji,Ej] ac- 
cording to 0{t) — 0{uj)e~'^'^* + C.C., with appropriate 
regularization implicit (Appendix O Sec. 



IIDl 



Having 

graphene in the Dirac cone approximation in mind, the 
latter equation contains the relevant degeneracies. The 
spin contribution as a degeneracy factor, gs, should be 
valid for typical magnetic fields (<15 T) when the Zee- 
man effect does not manifest. 

The first step is to project the Heisenberg EOM 
for the current onto the space of unperturbed single- 
particle states: we introduce the field operator ^cr(r, t) — 
X^ct ^cfO"(^)V'o[(r) (and the respective Hermitian con- 
jugate), where Ca.aicln a) is the annihilation (cre- 
ation) operator obeying fermionic anticommutation 

rules: {Ca.o-, cj^,,^,} = ^aa'^o-.cr' and {Ca^a} Cot' ,a'} = 

{^a cn cjj,/ ^/} = 0. The symbol a — (A, k, ...) specifies the 
single-particle state of the electron (or hole) and a = ±1 
is the spin variable. The kets \a.,a) = cj,, ^\0) represent 
eigenstates of Hq, and, therefore, the position representa- 
tion, {r\a,a) = "00:. ^(r), equals Eq. (jsj) at zero magnetic 
field or Eq. ^ in the presence of a transverse uniform 
magnetic field. 

The second-quantized form of the full Hamiltonian and 
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the current density operator is given by 
respectively, where 



3 = 



(8) 
(9) 

(10) 




is the cur rent d ensity of graphene in the continuum 
descriptionJiSIlIl We omit the spin dependence of the op- 
erators hereafter for clarity of exposition. 
We now define the generic operator, 



P«/3(<)^4(<)5/3W, 



whose EOM reads 
d 



dt 



^h^s P-is{t), Pocp{t) 

■7,5 



(11) 



(12) 



FIG. 2: Allowed interband transitions (vertical arrows) in 
graphene; a photon of energy hu produces an excitation from 
the lower to the upper Dirac, as long as hu > The tran- 
sitions conserve k and hence are said to be "vertical." For 
hid < 2/x, Pauli blocking forbids any (interband) transition. 
In practice, due to disorder (impurities, etc.), the interband 
conductivity can be non-zero even for fko < 2/i. 



With this notation, the current density along the x di- 
rection reads, 



where h^s = {llH S) are the matrix elements of the full 
Hamiltonian [Eq. (4|]. Solving for Pafi{t) gives directly 
the current density according to. 



Jx{t) 



~ E jlA',k,k(5l,k(i)5A',k(t)) 



(15) 



A.A'.k 



Ut) = Y,{oc\J^\|3)Pc.p{t), 



(13) 



The non-null matrix elements in Eq. ( 14 1 contributing to 



and hence the (yet non-regular) optical conductivity via 
Eq. ([7]). The regularization is the final step of the 
EOM approach needed for obtaining a fully-consistent 
conductivity (in particular, obeying Kramers-Kronig 
relations) .ESl The respective technical procedure is given 
in Appendix O 



the conductivity correspond to transitions between dif- 
ferent bands conserving the momentum k. These tran- 
sitions are said to be "vertical," and, in addition, since 
they connect states in different bands, they are refereed 
to as being interband- like (see Fig. [2]). 

Taking the dipole approximation, A{r) — > Ag, the 
Hamiltonian [Eq. (Isl] reads 



In the following section, we solve Eq. ( 12 1 explicitly in H = E] -^A(k)c^ j.CA,k 



+ C.C.) X 



the linear response regime (i.e., first order in the electric 
field) for any pair of quantum states q:,/3, in the absence 
of a magnetic field. The case of finite (nonzero) magnetic 
field intensity is left for Sec. |IID[ 



C. Graphene in a zero magnetic field 

The purpose of this section is to show the EOM method 
at work in the context of a simple problem, which allows 
us to derive well-known results. In the absence of mag- 
netic fields, the macroscopic electronic current follows the 
applied optical field, and thus only the longitudinal con- 
ductivity is nonzero. From symmetry considerations, we 
also have a^xi'^) = '^yyi^)- According to the statement 



A,k 



Jc.D.k.k^lk.o-^^.k.cr + (c 



(16) 



In the latter equation, ExCk.) = Xhvpk , and the sub- 
scripts c (v) denote electrons (holes). 

As described above, we need to compute the time evo- 
lution of the operator Pv,c,ki^) = k('')^c,k(0- Straight- 
forward algebra yields 



dt 



+evFAo (e + c.c.) ^.^ ["c(t) 



(17) 



Eq. ([13|, the relevant set of EOMs to be solved is de- where we have defined the occupation operator for elec- 
trons (holes) as ?t.c(ii)('^) = k(^)^c(t>),k(^)- A similar 

equation holds for Pc.v,k,tT which can be obtained by in- 
terchanging c i — > V. 



termined by the non-zero matrix elements of the current 
density. 



Defining (k, Ab,|k', A') = -{evF/A)jf 



A,A',k,k' 



and us- 



ing the wave functions Eq. (|3|, we easily find 



To proceed, we take the average of Eq. (17 1 with re 



JA.A',k,k' 



3k.k' 



(14) 



spect to the unperturbed Hamiltonian, Hq, and approxi- 
mate {ficit) — ny{t))o ~ {fic — ny)o. Both procedures are 
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consistent with an expansion of {t) up to first order in 
the parameter ^o- The solution of the above differential 
equation reads 



(18) 



with, 



and we have introduced an imaginary energy F by hand, 
so to account for disorder phenomenologically. The re- 
maining term Pv,c.Vi,a{^) can be obtained from the latter 
expression by making lo — w and F — > — F . From 
Eq. (15), the oscillator strength of the current density 



along the x direction Jx{oj) is seen to be given by 

(20) 

The longitudinal optical conductivity, Uxx, follows from 
Eq. 0, 



^intcr 



2, ,2 



e V 



+ (c ^ v) . 



(21) 



where npiE) ^ i/[(,iE-fj-)/kBT ^ -^j stands for the Fermi- 
Dirac distribution (/x is the chemical potential) . In deriv- 
ing this expression, we have used the relation Ex{uj) = 
iujAq. Taking the clean limit F — > and considering 
cj > and T = 0, one obtains the well- known result 



Re cr" 



vre 



' {hhj - 2\^l\) 



(22) 



The latter result is the T limit of Eq. 0. For 
photon energies higher than 2/^ (see Fig.|2|, the interband 
conductivity is essentially frequency independent (up to 
energies of ~2 eV) and equals 



7re 
'2h ' 



(23) 



which is nothing other than the universal conductivity of 
graphene mentioned in Sec. |l] For /i = 0, and contrary to 
ordinary semiconductors, there is no frequency threshold 
for interband transitions: according to Eq. (22 1, some 



interband transitions will always be available for a suffi- 
ciently high photon frequency. As a consequence, Drude's 
description will not suffice for a general description of the 
optical response of graphene. 

In addition to the interband transitions discussed here, 
there is an intraband contribution in graphene which can 
be appreciable for fi 0. This contribution comes from 
nonvertical processes (e.g., via collisions with phonons). 



not included in the Hamiltonian Eq. ([16|. This contri- 
bution gives the Drude response and reads^ 



Re (T 



intra 



2e' I 



F2 



(24) 



Interestingly enough, the latter result can be derived 
from a full quantum mechanical calculation by consider- 
ing a finite magnetic field intensity and taking the limit 
B — ?> in the end.^*^ This is because a magnetic field 
open gaps in the spectrum of a clean system, allowing 
for intraband transitions (see Sec. II D). A semiclassical 
calculation also leads to an equivalent result (Sec. II El. 



D. 



Optical conductivity of graphene in a magnetic 
field 



In what follows, we show that the EOM method can 
be employed to study the magneto-optical response of 
graphene along the same lines as in Sec. |II C[ The pres- 
ence of a transverse magnetic field in the Hamiltonian 
develops LLs, and hence we must start from the eigen- 
states given in Eq. (|6|. The latter defines the field oper- 
ator, ^'(r,t) = Y.,i,ky Cn,ky{t)^n,ky{r) (togcthcr with the 
respective Hermitian conjugate); the index n labels the 
degenerate LL with energy given by Eq. ([5|. The field 
operator can be written as 

^(r t) - — V ( hA-ii^) \ Jkyv^ 



(25) 



2L frl V *^^g^^('^)'^l"l(^) 

■^^E E (, I) 



This peculiar spinorial structure, with a single level being 
highlighted, is on the basis of non-standard features in 
the magneto-optical conductivity of graphene .-S'ii'^iiio, 



1. The longitudinal conductivity 

According to Eq. ([?]), the calculation of the longitu- 
dinal conductivity requires computation of the average 
value of the current density operator along the x direc- 
tion, 

= E E i^^h\M^',K)(^lkymn'Myit)) ■ (26) 



Using the LL wavefunctions [Eq. (|6|], we easily find the 
non-zero matrix elements to be, 



{0,ky\jx\±l,k'y) 

{n,ky\jx\n',ky) 



evp 



(27) 



i [sign(n')(5|„|_i,|„'| 



evp . 
' ~2A 

-sign(n)5|„|j„,|_i] 8ky,k'y , (28) 
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where in the last line n, n' ^ 0. These statements show 
that the optical transitions conserve ky and occur be- 
tween levels with indexes n and n! satisfying \n\ — \n!\ — 
±1. 

Two sets of transitions are thus allowed: intraband 
transitions, occurring within the same band, and, as 
in the absence of a magnetic field, transitions connect- 
ing LLs in the valence and conduction bands, which 
are interband- like. Transitions involving the zero-energy 
state n = can be considered either intraband- or 
interband- like, since the zero-energy state is shared be- 
tween electrons and holes. For the sake of simplicity in 
defining the set of EOMs, throughout, we classify transi- 
tions involving the zero-energy state as being interband. 

In order to clearly distinguish among the possible types 
of transitions, we define 



Cn for 71 > 

W|„l for n < , 



(29) 



ao 



for n = 



with the Hermitian conjugates following identical redefi- 
nitions. Note that with these definitions the subscript n 
in the operators take only positive integer values. 

a. Interband transitions — Using the field operator in 
the presence of a magnetic field [Eq. (25l[, and keeping 
track of just the interband terms for the moment, the full 
Hamiltonian takes the form 



H = ^ [Enclcn + E-nvlvn] 



evFA{t) 
evFA{t) 



^c\aQ + vlao + h.c. 



where A{t) = Ao(e~'"* +c.c.), and 



P 



(2) _ J 



(30) 



(31) 
(32) 



(Also, for clarity, we have omitted k y un der all the sum- 
mation signs.) The first line in Eq. (30) describes mass- 



less Dirac fermions in a transverse magnetic field and the 
remaining lines contain the electronic transitions among 
different LLs induced by the external electric field. 

The interband current density along the x direction 
can be recast into the form 



J: 



+ ^^^^ E ['Pn^ + 'Pn^ + li-^-) • (33) 

ri>l 

From the form of the current we see that there are two 
basic sets of EOMs to be solved: the first set refers to 



the time evolution of operators involving the zero-energy 
state (cjao, t'jao, and Hermitian conjugates), while the 
other set refers to higher energy LLs. Take, for instance, 
the operator p'^^ belonging to the latter set; as in the 



case of zero magnetic field (Sec. II C I, the commutator 
[iJ, P^^\ gives rise to (i) occupation number operators 
(v\j^^VnJ^\ and cJjC„), and (ii) a free evolution term, that 
is, the operator Pit^ itself. In addition, intraband terms 
with |n| - |n'| = ±2 show up, namely, cJjC„+2, i'l_iW„+i 
and q^qV^^ti,!- These terms do not originate real intraband 
transitions, since the respective current density matrix 
elements are null. 

We are now in the position to write the prototype 
EOMs governing the interaction of Landau quasiparti- 
cles with an external oscillating electric field. 



I dt 
h d 



1 



(34) 



^ ^^Pe = EiP, + -^evFA{t)[alaa ~ c\ci] , (35) 

where we have omitted the time dependence of the op- 
erators and defined Pc(v) — c{v)\ao. The remaining op- 
erators obey similar equations. [The EOM for P„ is 
obtained making P^^"* pl^-'and interchanging n with 



n + 1 on the right-hand side of Eq. ( 34 1 . As for P„, we let 
Pc Pv, El E_i, and Ci{c\) wi(w|) in Eq. 

To solve the above set of differential equations to first 
order in Aq, we proceed as in Sec. jll C| Taking the average 
value (...)o of each EOM with respect to the unperturbed 
Hamiltonian, Hq, the solution for each operator O can be 
written as {O{t))o = 0(tj)e-*"* -f- 0(-w)( 
oscillator strengths read 



where the 



= -^evFAn- ^. . ,T. ' (36) 



'-^'-(ri+l) — En — huj — iV ' 



(37) 



and where, as in Sec. |II C| we have added a imaginary 
energy F to account for level broadening. The solutions 
for P,\ '{uj) and P^,(w) can be obtained from the latter 
expressions as described below Eq. ( 35 1 . 



Combining these results and Eq. ( 33 I , we easily find 



J^(w) 



±evFY.lY.[^P,\H+^P^{u;)_ 

ky I n> 1 

-V2 ^Pciuj) + Pv{lu) + "c.c. term^ 



'} (38) 



where the summation over ky has been restored. This 
summation yields the degeneracy of the LLs J^k ~ 
A/{2Trlg). The last term in the above equation (i.e., the 
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E 



::t:z 

intraband 



■Nf + 1 
Nf 



-Nf 
-Nf - 1 
-Nf - 2 



6. Intraband transitions — The intraband interaction 
Hamiltonian reads 

Hl^^ = 'T^evpAit) J2 Rvn+i - 4c„+i - h.c] , (40) 



and the zero-energy operators {qq and Oq) are ab- 
sent given our classification of intraband transitions [see 
Eq. (29 1 and the following text]. The calculation follows 
identical steps to the interband conductivity and, hence, 
is not repeated. The final expression for the (regular) 
intraband diagonal conductivity reads. 



FIG. 3: Schematic of electronic transitions contributing to 
o'xxi^-^) of doped graphene in a magnetic field. In this exam- 
ple, Ef > El, and thus the last occupied LL, n = Nf > 1, 
belongs to the conduction band. Two types of transitions 
take place: (i) interband transitions, connecting LLs from the 
lower cone (valence band) with LLs in the upper cone (con- 
duction band), and (ii) intraband transitions within the up- 
per cone. Intraband transitions are limited to adjacent LLs: 
Nf — >■ Nf The figure shows the following interband tran- 
sitions: (a) the pair —Nf — >■ Nf -\- 1 and —Nf — 1 — >■ Nf, 
whose energy difference is -Ejv+i -I- -Ejv (the lowest interband 
energy; note that transitions — A'^f — 1 — >■ Nf are forbidden be- 
cause n = Nf is occupied); and (b) the pair —Nf — 1 — >■ Nf+2 
and —Nf — 2 — ^ Nf -\- 1. The respective energy difference is 
En+i+En+2 (the second lowest interband energy difference), 
and in this case both transitions take place. Transitions with 
higher energy differences are not represented. 



c.c. term) is obtained taking the complex conjugate and 
making uj — > —uj of all the previous terms. 

The final expression for the longitudinal (interband) 
conductivity is derived in two steps: (i) dividing the 
Eq. (38 1 by Exi^^) [Eq. (m)[, and (ii) undertaking ap- 



propriate regularization to remove the divergent factor 
l/uj, 



N^ 



-z^(H-,5„,o) J2 



n=0 



a=±l 



_£;_(„+!) - En - En - a{hu} + iT) 

+ {n o n + 1)] . (39) 

The above expression is analytic in the upper-half plane 
and finite at lu — 0, thus obeying Kramers-Kronig causal- 
ity relations. (We refer to Appendix [C] for the derivation 
and physical grounds of the regularization procedure.) 
Note that, as usual when dealing with low-energy theo- 
ries, a cutoff energy i?cut of the order of the bandwidth 
must be considered for consistency; we take n < N^, 
with Nc — int[(i?cut/-6'i)^], where int[...] denotes the in- 
teger part. Nc varies roughly as 10'^ with B in teslas. 
Within the physical relevant range for Ecut, these sum- 
mations converge quite rapidly; the figures in the present 
work have Ecut « i — 2.7 eV. 



^intra 



(^) 



« Q = ±l n=l 



En+1 — En 



^ nF[En+l] - llFjEn] 

En+1 - En - a [huj + iV) 

+ {En — i> —En A En+1 —En+l)] 



(41) 



The full longitudinal conductivity a^xi^) is given by 
adding its interband and intraband counterparts, that is. 



Eqs. (39 1 and (41 1, respectively; straightforward algebra 
yields 



(^xxil^) 



with E„ 



e 
~h 



Na 

E 



A XX 

nrn 



np^En) - npiEn 



-Na 



iE„ 



hjjj + En 



iT 



(42) 



Er, — E.-n . and where we have defined the 



longitudinal matrix elements 



A XX 

nm 



;2 



-(1 -I- 5ni,0 + 5n,o)5\ 



i»i,±i' 



(43) 



Equation (42 1 is the main result of the present sec- 
tion. It coincides with Eq. (7) in Rcf. Jjj obtained via 
a Green's function calculation in the bubble approxima- 
tion and, also, with a Kubo formula calculation within 
the Dirac cone approximation (see Appendix|C|. We note 
in passing that, on top of the interband and intraband 
contributions discussed here, there is a correction arising 
from phonon-electron coupling. At low temperatures and 
zero field, this correction is expected to be small."^ At a 
high magnetic field, though, a recent calculation shows 
that phonon energy peaks split the LLs nearbyj32l which 
can lead to a measurable signature in magneto-optical 
experiments. 



2. The general properties of axx{^) 

In what follows, we overview the main features of 
graphene's longitudinal magneto-optical conductivity, an 
essential step to understanding the Faraday rotation in 
graphene (Sec. III). 

a. Low electronic density — At a low electronic density, 
more precisely, for \Ep\ < Ei, no intraband transitions 
can take place. Because the LL energy scale in graphene 
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FIG. 4: Longitudinal magneto-optical conductivity as a func- 
tion of the photon energy for a field of 7 T, zero chemical 
potential, T = 17 K, and F = 6.8 meV (~ 79 K). The hori- 
zontal dashed-dot (black) line marks the graphene's universal 
ac-conductivity background [Eq. |TJ]. 



is relatively high (e.g., Ei ~ 36 meV for a field of 1 T), the 
magneto-optical conductivity is fully driven by interband 
transitions even close to room temperature. 

Figure |4] shows a plot of Eq. (42 I for zero Fermi energy 
and a magnetic field of 7 T: a sequence of absorption 
peaks, corresponding to the maximum of the real part of 
each term in Eq. (39), fiLj ~ Ei, £'2 — E3 — E-2, etc.. 



is clearly observedjsee Eq. ( 49 1 and text thereafter] . The 
conductivity never vanishes, even though the concentra- 
tion of carriers is low (Ep — 0), a genuine signature of 
graphene's LL structure.'^ 

The contributions from different interband transitions 



[Eq. (39 1] partially overlap at a high frequency, with the 
effect that the real part of axx{^) displays the so-called 
Shubnikov-de Haas oscillations around the universal ac 
optical conductivity of graphene, CTg (the imaginary part, 
in turn, oscillates around 0).'Sfi^The semiclassical con- 
ductivity is null, on the other hand, thus failing to de- 
scribe the magneto-transport in neutral graphene. 

b. High electronic density — Away from charge neutral- 
ity, more precisely, for \Ep\ > Ei, the picture is more in- 
volved; intraband transitions can now occur, while some 
interband transitions will be blocked. We take T — 
and, without loss of generality, assume that Ep > 
(similar conclusions hold for holes); direct inspection of 
Eq. ( [4l| shows that a single type of intraband transition 
is allowed, whose contribution to the optical conductivity 
reads 



2ihvp 



huj + iV 



In the above formula. 



^Awat = En j^i — E^ 



(44) 



(45) 



denotes the intraband gap, with Np being the index for 
the last occupied LL. 



Let us first consider the limiting case when the energy 
gap /^Nf is larger than the level broadening, A^r^ > T. 
The latter typically happens at high magnetic fields and 
not too high Fermi energies; in this limit, the real part 
of Eq. ( 44 1 displays a maximum at w ~ Awatj^ , with an 



intensity falling off as B/Aujn^ 



2eBvj. 



nTAu 



Nf 



X cr„ 



(46) 



The intraband magneto-peak, Eq. (46), is the lowest fre- 
quency peak in the absorption spectrum of graphene 
with Ep > El] its magnitude increases with increas- 
ing Fermi energy and/or magnetic field intensity. An 
example of an intraband absorption line occurring at 
uj ~ AuJNp is shown in Fig. [Sj In that case, the param- 
eters correspond to Ajv^ — 22.6 meV and F = 6.8 meV, 
and hence A^r^ > F. Some points are worth men- 
tion: (i) the intraband contribution to the conductivity 
[Eq. (44l[ dominates at low photon frequencies; and (ii) 



the curve for Re axx{^) shows that the remaining absorp- 
tion peaks are found in the higher frequency part of the 
spectrum, above the threshold for interband transitions, 
hu > E^F + Enf+1- (Note that, at a low magnetic field 
and/or high Fermi energy, the level spacing between ad- 
jacent LLs is so reduced that iJ^Vp ~ E^j^j^i — Ep, and 
thus one recovers the condition found earlier, namely, 
huj > 2Ep.) Such interband peaks cause Shubnikov-de 
Haas oscillations despite the finite electronic density. 

For a general relation between the broadening and the 
energy gap A , the maximum for the intraband peak 
occurs at 



(47) 

When Aojnf — r/(-\/3/i) (typically the case for a very 
high Fermi energy and/or low magnetic field), the intra- 
band conductivity is maximal at null frequency, with an 
intensity given by Eq. (46) multiplied by a factor of 2. 

The regime A^v^ ^ F is illustrated in the bottom panel 
in Fig. [T2j Two magnetic fields are considered, at a 
fixed Fermi energy, Ep = 0.3 eV, with Re axxi^) be- 
ing represented by the solid lines. When B = 7T (left- 
hand panel), although a considerable number of levels 
are occupied {Np = 9), one has A^^ ~ 1.4r, which, ac- 
cording to Eq. (47), corresponds to a maximum of the 
longitudinal conductivity at a; ~ Auj^f- This is in- 
deed confirmed by the numerical calculation shown there. 
Decreasing the magnetic field down to B —3 T (right- 
hand panel), reduces Ajv^ (recall that the LL energy 



varies as ^ y'B), which in turn increases the num- 
ber of occupied levels to Np = 22. As a consequence, 
AjVj? — 0.67F, and the maximum of the intraband peak 
is seen to be shifted to zero frequency, again in accor- 
dance with Eq. (|47|). 



Given the intrinsic large cyclotron gap of graphene. 
El, the intraband contribution [Eq. (44l[ controls the 
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FIG. 5: The longitudinal conductivity as function of the pho- 
ton energy for Ef = 0.2 eV. Other parameters as in Fig. [4] 
The solid horizontal (black) line shows graphene's universal 
ac-conductivity background [Eq. ([TJ]. 



magneto-optical response of this material in the mi- 
crowave region up to terahertz (THz) frequencies in sam- 
ples with a finite electronic density {Ep > Ei). 

The interband contribution, on the other hand, is im- 
portant both in samples with a low electronic density, 
Ep < El, where it determines the full magneto-optical 
response (discarding the effect of phonons as discussed 
above) , and in samples with arbitrary carriers concentra- 
tions, for photon energies above the threshold for inter- 
band transitions, hu: — Eiq^ ~\- E^^^i (typically within 
the near-infrared region). 

The positions of each interband peak can be obtained 
from Eq. (47), with Awn^ replaced by 



{En+i+En)/h, 



(48) 



with the constraint n > Np- At finite electronic densi- 
ties, Np > 1, typically one has HAiln ^ T, and thus we 
arrive at the following useful approximation 



inter (n) ^ a q 



,n> Np. 



(49) 



For not too small fields, B > 0.1 T, the cyclotron gap 
£;i ~ 36 X ^ meV-T-i/2 jg ig^j-g^j- than the LL broad- 
ening, and thus, in practice, the latter statement can be 
generalized to include the case of Np = 0. 

For general parameters, the intensity of each interband 
peak is no longer given by a simple expression, because 
many interband transitions can contribute to the spectral 
weight close to each of the resonances uj ~ AJ7„. As 
a result, as uj varies, the real part of axxii-^) oscillates 
around a constant value of about cTg. Examples are shown 
in Fig. |4] for Ep ^ and in Fig. |] for Ep = 0.2 eV. In 
the first case, we have Np = and therefore all the 
observed peaks are interband-like. The second case has 
Np = 4 and therefore one intraband peak is observed, 
corresponding to transitions n = 4— )-n = 5, at low 



photon energy, whereas the interband peaks appear at 
energies hw > 2Ep = 0.4 eV. 

We finally remark that, as long as not too low magnetic 
fields are considered {B < 0.1 T), the above considera- 
tions are valid even close to room temperature (e.g., for 
_B = 1 T, the first LL corresponds to a thermal energy of 
420 K). 



3. The Hall conductivity 

The Hall optical conductivity of graphene, axy{uj), fol- 
lows directly from Eq. choosing i — y, j — x, we 
obtain 



Jy(uj) 

-9s9v X -fj- 
Ex{uj) 



(50) 



where we have invoked graphene's sixfold crystallo- 
graphic symmetry to write axyi(^) = —ayx{uj). The cen- 
tral quantity to be computed this time is the average 
value of the current density operator along the y direc- 
tion; using Eqs. ([6| and ( p5| , we get 



Jy{t) = -evp^{n,ky\jy\n' ,ky){c\,^^{t)cn> ,k-^{t)) 



The non-zero matrix elements read 

. evp 



{Q,ky\jy\ ± l,ky) 



evp 



(51) 
(52) 



{n, ky\jy\n' , ky) = - ^ [sigU ( Tl' ) <5| « h 



2A L"'''"^"' 

+sign(n)(5|„|j„,|_i] , (53) 

(plus respective complex conjugates) where, in the last 
line, n,n' ^ 0. Omitting the summation over ky, the 
total current density reads 

Jyit) = '^T^^'^^F + v\aQ - h.c.^ 



2A 
1 

2A 



evp 



h.C. 



(54) 



TI>1 



The EOMs resemble those derived for the longitudinal 



conductivity [Eqs. (34)-(35l], the reason being that the 



current matrix elements in the x and y directions are the 
same except for phase factors [compare Eqs. ( |27[ ) and 
([28]) with Eqs. ^ and ([53])[. The final formula (after 
regularization) yields. 



rog/ Kmi np{En) - np{Em) 

"xy \^) — ~r /, ^ T- r~B I -p ' 

n^rn——Iyc 



with matrix elements related to AJj^ [Eq. ( [42| )[ ac- 
cording to, 

A^m = iA^^(5|„|j„|_i - 5\rn\-l,\n\) ■ (56) 
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the interband gap, hAft^p. The only exception is in- 
deed the interband transition —Np — >■ Np + 1 because, 
contrary to interband transitions involving larger energy 
differences, it cannot be canceled by the other member 
of the pair, n — —Np — 1 A m = Np, since the latter is 
forbidden via Pauli blockade; a schematic picture is given 
in Fig. [6] 

The extremum points of the real part of the Hall con- 
ductivity occurs at zero frequency, w = 0, and 



(59) 
(60) 



FIG. 6: Schematic of electronic transitions contributing to 
the Hall conductivity of doped graphene in a magnetic field. 
Contrary to the longitudinal conductivity (Fig.[3|, symmetry 
implies that only interband transitions involving the smallest 
energy diflference, hAQ^p = E^p -\- iSjVj^+i, contribute to 
axy The remaining interband transitions (Af2„, with n > 
Nf) come in pairs whose contribution to the Hall current 
mutually cancel as explained in the text: an example of a 
pair of interband transitions that cancel is shown in zig-zag 
arrows. Note: The schematic picture is strictly adequate for 
Np > 1; the case of Np = admits a single type of electronic 
transition, namely, n = — >■ n = 1. 



where we have considered T/h < Aujj\jp [see Eq. (46 1 and 
text therein] and made use of F/fi <C Aftpfp. The latter 
consideration is true in virtually all situations except for 
graphene at a low electronic density and small magnetic 
field B. Within the same accuracy, the Hall conductivity 
at w = reads 



Re axy{0) 



1-5 



Nf,0 



1 + o\ f ^eBv 



whereas at the point lo™^''''^ it is given by 



XCTg, 

(61) 



Likewise Uxxi^)- the result for the Hall conductivity 
based on the EOM method coincides with the result ob- 
tained using Green functions calculations .^D 

Symmetry considerations imply that only two terms 
contribute in general for the zero-temperature Hall con- 
ductivity, and hence the formula Eq. ( 55 1 can be consid- 
erably simplified. The first term is the intraband contri- 
bution and reads. 



' xy 



1 



h 1% {hiu + ivy - n?Auj% 



(57) 



and the second is interband-like, connecting electronic 
states with n — —Np and n — Np + 1, and reads. 



1 + 5np,o 



(58) 



A single interband transition play a role in setting the 
Hall conductivity, even for zero Fermi energy. This is 
at odds with the situation for axx{^), where many non 
equivalent interband transitions contribute to the optical 
spectral weight. To understand this peculiar feature of 
axy{uj), let us consider the second lowest interband res- 
onant energy, namely, AE2 = Emp+2 — E-iqp-i: there 
are two distinct sorts of interband transitions n ^ m 
involving such energy difference, namely, the pair ni — 
-Np-2 A TOi = Np+1 and = -Np-1 A ma = Np+2, 
whose Hall matrix elements read, A.^\rni 

\xy ^ 



and 



^nams ' respectively. When substituting into 
Eq. (|55|, these contributions cancel each other at T = 



because A^^„^ 



a; 



The same argument applies to 
all transitions involving an energy difference larger than 



Reaxy{u:'l''^)c^F, 



± /I - 6Npfi\ ( eBvj, 



Auj 



and for it reads 



Reaxyioj'T"'] 



F. 



An. 



1 



ttF 



eBvl 



X cr„ 



(62) 



X (T„ 



AQnp J \ ttT 

(63) 

where we have defined = ±huj/{fiuj ± F). The in- 
tensity of the Hall peaks dependence on the magnetic 
field intensity B is the same as for the longitudinal (in- 
traband) peaks [Eq. (46l[, i.e., as ~ Vi?. Also, similarly 
to <Txx{'-^)i in doped graphene with Np > 1, the inter- 
band peak is very low compared with the intraband Hall 
peak for AwjVp ^ Ail^p. We, finally, remark that the 
anomaly associated with the zero energy LL is present in 
all the latter expressions via the factor 1 + Sj^p Q. 

Figure [7] shows the Hall conductivity of graphene at 
a high magnetic field {B — 7 T) for Np — (top) and 
Np — 4 (bottom), corresponding to neutral and highly 
doped graphene samples, respectively. The main charac- 
teristics of Reaxy{u}) can be explained using Eqs. (61)- 
(63 1. In particular, for doped graphene, the spectral 



weight concentrates around two well-separated parts of 
the spectrum: (i) an intraband-dominated region (n = 
4 — > n = 5 ), at low photon energies, with a maxi- 
mum (minimum) intensity occurring at hu:^ ~ A4 -f F ~ 
30 meV (/iw_ ~ A4 — F ~ 16 meV ) [intensity equal 
to ~ lOe^/h (~ —20e^/h), in accordance with Eq. (62l[, 
and (ii) an interband-dominated region n = — 4— >n = 5, 
at high photon energies, with a maximum (minimum) 
intensity occurring at ~ hAD,4^ -t- F ~ 413 meV 
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FIG. 7: Hall conductivity as a function of photon energy for 
Ef = (top) and Ef = 0.2 eV (bottom). In both plots 
T = (other parameters as in Fig.g. At zero Fermi energy 



(top), 



originates in a single type of interband transi- 



tion, centered at hui ~ Ei 96 meV, and ther efor e cannot 
be described by a semi-classical treatment [Eq. (58l]. When 



Ef = 0.2 eV (bottom), the first four LLs are fulfilled, which 
results in a classical intraband contribution [Eq. ( |57[ l], cen- 
tered aX huj ~ E5 — E4 ~ 23 meV, and a single interband tran- 
sition [Eq. ([58|[ centered at huj ~ E5 + Ei ^ 0.4 eV ~ 2Ef. 
The latter is shown in the inset. 



{hfl^ ~ hAfti - r ~ 400 meV) [intensity equal to 
~ O.SleV^i -OMe^/h), in accordance with Eq. 

Dependence on the Fermi energy — The variation of 
conductivity with the Fermi energy reveals other peculiar 
feature of 2D systems: Hall quantization.l^SlMl Figure [s] 
shows the formation of plateau in the static (or dc) Hall 
conductivity, a^yiO), a direct evidence for discrete energy 
levels. In conventional 2D electron gases, the widths of 
such plateau are constant (the LLs energy scales as n), 
whereas in graphene the plateau's width decreases with 
increasing Fermi energy (the LLs energy scales as y^n). 
As for the steps heights, they are equidistant in graphene, 
A(Txy{0) ~ Ae'^/h, even when crossing Ep = 0, whereas 
in conventional 2D systems the step from the first elec- 
tron LL {n ~ 1) and the first hole LL (n = — 1) is twice 
the value of the remaining steps (a manifestation of the 
zero-energy LL graphene anomaly). 

The Hall conductivity quantization rule for graphene 
can be readily obtained by adding the intraband and in- 
terband Hall conductivities, 



4e2 / 1 



(64) 



where we have used T ^ Ei in order to simplify the 
denominators of Eqs. (57) and (58). Despite the filling 
factor, V — ANp + 2, being an integer number, there 
is no complete correspondence with the conventional 2D 
IQHE, for which a^y = —'^e^Np/h; an extra 1/2 factor 
due to the contribution of the zero-energy state, shared 



FIG. 8: The dc Hall conductivity as a function of the Fermi 
energy. The parameters are T = 17 K and F = 0.68 meV. 
The plateaux show Hall quantization values according to the 
theoretical prediction for massless Dirac fermions [Eq. (64|]. 



by both electrons and holes, shows up, which must be 
taken separately, making v always even — this is known 
as the anomalous IQHE and is a hallmark of chiral 
massless fer mion s. The anomalous IQHE was predicted 
theoretic all jEEHl and measurecP^^ in the early days of 
graphene. 



E. The semiclassical solution 

Within the semiclassical approach, the spinorial na- 
ture of the electrons' wave function is immaterial. On 
the other hand, the massless nature of the spectrum in- 
validates a straightforward Drude-like approach^"^ to 
the calculation of the transport coefficients, and Boltz- 
mann transport theory is required, since in this formula- 
tion the central quantity to be computed is the deviation 
of the momentum distribution function from the equilib- 
rium Fermi distribution. 

In the semiclassical regime (that is, a high electronic 
density and/or low magnetic fields), the physics of the 
Hall effect can be explained in terms of Boltzmann's the- 
ory of transport, where the electric current is given, in 
the case of graphene, by: 



J = 



(65) 



with spin and valley degeneracies included, and where 
.g(i?, k, w) = (7k, is the deviation of the carriers' (elec- 
trons or holes) distribution function from the equilib- 
rium Fermi distribution, /o(e), e is the charge of the 
carrier, the static magnetic field B is considered to be 
perpendicular to graphene's surface, cj is the frequency of 
the electromagnetic field, and the carrier's velocity reads 
■Wk = {vx,Vy) = vp{cos6,s'm6). In the presence of both 



12 



an electric and a magnetic field, the distribution gt is the 
solution of the equation'*'^ 

- eE • = \- — + -{vuy.B) ■ Vkffk , (66) 

oe Tk OT ft 

where we have employed the the standard relaxation 



approximation,^ i.e 



dt 



9k 

Tk 



(67) 



where Tk is the relaxation scattering time, E — 
{Eq^x, Eo^y) is the electric field, and Vk is the gradient 
operator with respect to the momentum k. Writing (7k 
as. 



B=7T,Ep = 0.1 eV(Np=1) B=7 T,Ep = 0.15 eV (Np=2) 




100 200 300 100 200 300 
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5k 



*k - A, 



(68) 



and noting that (t^k x B) • Vk5k — fk • (B x Vk5k), 
Eq. ( 66 1 can be solved exactly, where the vector A needs 



to be determined. Solving Eq. (66 1, the components of 
the vector A = {Ax,Ay) are obtained in the form 



A.r. 



Ar,, 



(1 - iujTk)Ex - T]f,UJcEy 
(1 - iujTk)Ey + T]^UJcEx 



where 



(l-ZWTk)2+w2r2 



uJc^ev%B/\EF\ 



(69) 
(70) 

(71) 



is the graphene's cyclotron frequency, and E^^y) is defined 
as 



Ex(y) - '-eEo^xiy)Vxiy)-Q^ 



(72) 



Introducing gk in Eq. (65 1, and assuming T = 0, we 



obtain the components of the conductivity tensor, which 
read 



e^2\EF\Tk, 



1 - iujTkf 



h h {1 



(73) 
(74) 



Note that setting Wc = in Eq. ( 73 1 leads to the semi 



classical longitudinal conductivity at zero field mentioned 
in SecHTCl 

Validity of the semiclassical calculation — The results 
presented so far demonstrate the reliability of the Boltz- 
mann approach in regions of the spectrum where the op- 
tical weight is mostly due to intraband transitions. This 
is borne out in Fig. [s] [Fig. [t] (bottom)], where axxi^^) 
[axyito)] is plotted as a function of hoj, for i? = 7 T and 
Ep = 0.2 eV: the agreement between the real part (imag- 
inary part) of the quantum calculation shown by the 
solid (blue) line [dashed-double-dotted (red) line] and the 



FIG. 9: The real part of the longitudinal conductivity is plot- 
ted as a function of the photon energy for Ep = 0.1 eV (left) 
and Ef = 0.15 eV (right). In these plots, B = 7 T and 
r = 6.8 meV. The dashed (red) represents the semiclassi- 
cal result [Eq. ( |73[ )[ and the s olid (blue) line represents the 
EOM quantum solution [Eq. (42l]. Note that, in the right 
panel, there is no interband peak n = — 1— )-m = 2, at 
Ef ~ 230 meV, and the peak at htu ~ 300 meV loses half of 
its intensity because the n — —3 — >■ m = 2 transitions get 
blocked when the Fermi energy crosses the LL with n — 2. 



semi-classical calculation shown in dashed (green) curve 
[dashed-dotted (orange) line[ in these figures is confined 
to energies huj < 2Ef- For high photon frequencies — 
more precisely, above the threshold for interband transi- 
tions, huj ~ 2Ef — the conductivity cannot be described 
by Boltzmann's transport theory. 

The fine agreement observed at low photon energies 
is not accidental and ceases to occur only for a very 
low Fermi energy. To see why, we note that Eqs. (44 1 



and (|57| (intraband conductivity) and Eqs. (73 1 and (74 1 
(semiclassical conductivity) coincide, upon identification 
of the intraband energy gap A , with the cyclotron en- 
ergy hujc- This identification is justified when a sufficient 
number of LLs are filled. In fact, expressing the Fermi 
energy as Ef = {Hvf /Ib)V2N* , we obtain A^r^ — )■ hojc 
provided that 



(75) 



Noting that Nf = int[iV*], we then see that the latter 
limit is achieved when N* ^ 1, as anticipated. 

For the parameters in Fig. |5] (see also Fig. [7j bot- 
tom), even though only a few LLs are fulfilled, i.e., 
Nf — 4, the values of A^r^ and hujc are quite similar, 
Ajv^ = 0.0226 eV and fiWc = 0.0230 eV, explaining the 
consistence between the two theories in describing the in- 
traband electronic transport. In practice, only for a very 
low Fermi energy and/or a very high magnetic field, such 
that Nf — 0, does the semi-classical calculation fail to 
describe the conductivity in the whole optical spectrum. 
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since all transitions are interband-like in this case. Re- 
markably, already for a single occupied LL, Np = 1, the 
semiclassical calculation provides a reasonable descrip- 
tion of the optical conductivity, as long as one remains 
inside the portion of the spectrum where the interband 
processes have little or no weight, that is, hw < Ei + E2 
(see Fig. [9]). We note again, however, that the intraband 
region extends for a large range of frequencies given the 
large intrinsic cyclotron gap of graphene. 

In summary, the validity of the semiclassical calcula- 
tion is bound to photon energies below the interband 
threshold, fiuj < E^f^ + E^p+i, and for not a too low 
Fermi energy, Np > 1. For the parameters used in 



Figs. |5] |7] (bottom), |9] and [TT] we list in Table I the 
corresponding values of Apfp/hiOc- These figures have 
Np > 1 and hence the semi-classical conductivity agrees 
well within the far-infrared part of the spectrum. For 
completeness. Fig. 9] shows the real part of axxi^^j) for 
Np ^ 1 (left) and Np = 2 (right). The former has 
An^/HuJc — 0.86 and hence the semi-classical calcula- 
tion is only partially accurate. In particular, it under- 
estimates the maximum intensity for intraband light ab- 
sorption. The right panel in Fig. |9] with Np = 2, has 
Amp/IvjJc — 0.99, which explains the excellent agreement 
between the two curves in the intraband region, huj < 0.3 
eV. 



B (T) 


Ef (eV) 


Nf 




1 


0.30 


68 


0.9990 


2 


0.30 


34 


0.9954 


5 


0.30 


13 


1.0066 


7 


0.20 


4 


0.9837 


7 


0.15 


2 


0.9933 


7 


0.10 


1 


0.8629 



TABLE I: Values of several relevant quantities related to the 
numerical simulations given in Figs. [517191 , and[TT| The agree- 
ment between the semiclassical calculation and the quantum 
intraband expression comes from the similarity between Ajv^ 
and fwjc- 

Having presented the calculation method for the 
magneto-optical properties of graphene based on the 
EOM method, we we now turn to study of the Faraday 
effect. 



ellipticity, and Faraday rotation angle. These quantities 
depend on the frequency of the impinging light, w, mag- 
nitude of the (transverse) magnetic field, B, scattering 
mechanisms (i.e., level broadening, T), temperature T, 
and Fermi energy, Ep, via the magneto-optical conduc- 
tivity tensor of graphene derived in Sec. [ll] 

Our theoretical results are tested against experimental 
data measured recently by Crassee et al. using graphene 
samples with a high electronic density .12^ The limit of 
a low electronic density is studied in Sec. |III C| where 
the Faraday rotation angle is shown to display quantum 
jumps as a function of the Fermi energy. 

Finally, in Sec. [Ill D| an experimental setup is proposed 
that is able to greatly enhance the Faraday rotation angle 
in the entire optical spectrum. 



A. Faraday rotation in graphene 

We now solve the problem posed in Fig. |10| considering 
only a single graphene sheet separating two dielectrics. In 
what follows, we assume that graphene is deposited on 
top of a lossless dielectric medium (i.e., fully transpar- 
ent to impinging light), of relative permittivity e^. The 
generalization of the problem to the case of a lossy dielec- 
tric poses no difficulties, except for the introduction of a 
complex index of refraction associated with the dielectric 
medium. We further assume that the incoming electro- 
magnetic field is linearly polarized along the x axis and 
propagates along the z direction, as shown in the diagram 
in Fig. [lOl that is. 



= e^Eie'^"'"'^*'^ 



(76) 



such that q — ^JT^u^j c. 

Due to the optical Faraday rotation of the plane of 
polarization of the electric field, both the refiected, W , 
and the transmitted, . fields acquire a finite component 
along the y-direction; that is. 



(77) 
(78) 



where k — lo /c. For this problem. Maxwell's equation for 
the electric field reads (in MKS units) 



III. THE FARADAY EFFECT IN GRAPHENE 



iujliQ5{z) (JijEj + Lu^erfiEi = , (79) 

j=x,y 



We discuss the transmission of electromagnetic radia- 
tion between two dielectric media separated by graphene. 
The scattering geometry is given in Fig. |10| where the 
transverse magnetic mode is chosen as a particular ex- 
ample. Since we are interested in a normal incidence, 
there is no distinction between the transverse magnetic 
and the transverse electric modes. 

The present section is organized as follows: in 
Sec. [Ill A] we derive general expressions for transmission. 



where Ei is the i-component of the electric field (we have 
i — x,y), iiQ is the vacuum permeability, and aij are the 
components of the magneto-optical tensor of graphene 
(see Sec. [ll]). The boundary conditions at the substrate- 
graphene-air interface are the continuity of the tangential 
components of the electric field at the surface of graphene 
(z = 0), 



(80) 
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graphene 



FIG. 10: Schematic of the Faraday effect: an electromag- 
netic wave polarized in the xy plane (transverse magnetic 
mode) and traveling in the positive 2 direction passes through 
a graphene film subjected to a transverse magnetic field B. 
In this case, graphene is adhered to a substrate (typically 
8102), but the experiment can also be made with suspended 
graphene. The transmitted field sees its plane of polariza- 
tion rotated by an angle Of and acquires a certain degree of 
ellipticity. 



and (the derivatives are evaluated at z = 0) 



dEf 



dE] dET 



dz dz 



dz 



(81) 



wherethe last condition was derived from integrating 
Eq. ((79| in the interval z e [0^,0+] and I = x,y. Cal- 



culation of the transmitted intensities becomes easier to 
perform if we rewrite the boundary conditions in terms 
of circularly polarized waves: 



-2qEl + (fc + q)E^^ = -^nwa^E] 



(82) 



where E± = E^^iEy and cf± 
resentation the two circular polarizations decouple from 



axx^i(^xy- for in this rep- 



each other. From Eq. (82 1 the transmission amplitudes 
follow in the form: 



*^=E^ 



1 + ^/£r + CMoCTzp 
The transmittance can be written as 
1 



\t±\e 



ie± 



T{B) = 



{\U? + \t-\'). 



(83) 



(84) 



where the factor 1/2 comes from the proper normaliza- 
tion of circularly polarized waves (omitted in the defini- 
tion above, for simpUcity of writing) and the factor 1 / 
is due to flux conservation. Faraday's rotation angle Op 
and the ellipticity 5 are given by**^ ^-^ 



s = 



(85) 
(86) 



respectively. From Eq. (83 1, 9p is given in terms of the 
conductivity a±, since 



arctan 



1 + -/fV + C/ifT^ ' 



(87) 



where (J± = (7'± + i(J± , and (Tj_ and cr^ are the real and 
imaginary parts of a± , respectively. Explicitly, we have 



O'i = (c^L T (^xy) + i {^xx ± 



xy) 



from which follows the approximate expression 



1 



' xy : 



(89) 



where we have assumed that 6*^- < 1 and that 1 -I- y/e^ ^ 
cfiQa':^. The latter assumption is the more stringent of 
the two. For comparison, in the numerical studies shown 
in Fig. |11[ we represent both the exact and the approxi- 
mate results for Op, S, and T. This allows us to check the 
validity of the approximate r esul ts. Discarding terms of 
the order of (c/ioCip)^ in Eq. (84 1, we obtain an approx- 
imate expression for the total transmitted hght in the 
form 



T{B) 



4^/1. 



(1 + V^rY 



1 



2c^o 



1 + 



(90) 



Within the same degree of approximation used to derive 
Eq. ( 89 ) , the ellipticity is given by 



2c^ 



'xy ■ 



(91) 



The validity of these approximations depends on the 
photon frequency, as can be seen in Fig. [TT] In what fol- 
lows, the exact expression is used in all numerical studies. 

In our simulations of the Faraday effect, we assume 
broadenings of the order of 10 meV. Our assumption is 
consistent with the values found in pump-probe exper- 
iments performed in exitaxial and exfoliated graphene 
samples''^ and in infrared spectroscopy studies of the 
Drude conductivity of graphene.^^ 



B. Fit to experimental data in the high-density 
regime 

Figure [T2| shows fits for two sets of experimental 
data for 9p^^ measured when electromagnetic radiation 
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FIG. 11: Faraday rotation angle (in degrees), normalized 
transmittance, and ellipticity of electromagnetic radiation 
passing through graphene subjected to a perpendicular mag- 
netic field. The graphene sample is assumed to have a finite 
electronic density, Ep — 0.3 eV, and to be on top of Si02 
(er = 3.9). Top six panels: Simulation of Op, T{B)/T{0), and 
5, considering a broadening of F = 7 meV. Bottom six panels: 
Simulation of the same quantities as above for F — 3.7 meV. 
In all panels, dashed lines correspond to approximate calcu- 
lations, as given by Eqs. and (|9T1), and T = 17 K. 
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FIG. 12: Faraday effect in doped graphene. Top: The Fara- 
day rotation angle (in degrees) when graphene is grown on 
silicon carbide. Fit to the experimental values oi Of, at a 
magnetic field of B = 7 T (left) and B = 3 T (right), us- 
ing the semiclassical approach [dashed (green) line] and the 
full quantum calculation [solid (red) line]. Parameters are 
Ef = 0.3 eV, r = 10.5 meV, T = 6 K, and = 4.4. Bot- 
tom: Theoretical optical conductivity ]Eqs. (42 1 and (|55[l] 



for the same parameters used to fit the experimental data: 
B = 7 T (left) and B = 3 T (right). 



passes through graphene epitaxially grown on siHcon car- 
bide (data taken at a temperature of 6 K). According 
to the experiments by Crassee et al.]^ it was possible 
to produce a single graphene sheet grown on the Si- 
terminated surface of 6H-SiC (the sample underwent H 
passivation of the Si dangling bonds, resulting in quasi- 
free-standing single-layer graphene) . Two sets of experi- 
mental data are shown in Fig. 12 (top), corresponding to 
two magnetic field intensities, i? = 7 T and i? = 3 T. 

In the course of the experiments it was found that 
the bare substrate did not reveal any Faraday effect, 
and therefore the measured rotation angle is intrinsic 
to graphene. This statement is confirmed by the model 
developed in Sec. |II1 A[ ARPES measurements on the 
sample used indicated a Fermi energy of the order of 



Ef ~ 0.34 ± 0.01 eV. 

In order to fit the data we have used Ep = 0.3 eV. We 
do not expect a perfect fit because we are considering a 
lossless dielectric. Nevertheless, the fit is fairly accurate, 
given the simplicity of the model. Moreover, the value of 

was set to 4.4 which is not the relative permittivity of 
SiC and must be understood as an effective number, given 
that the experimental data were taken with epitaxially 
grown graphene. Although the calculation in Sec. |II D| 
does not include this fact explicitly, the fits are satisfac- 
tory, for they reproduce the main features of the exper- 
imental data: a decrease in 9p with the photon energy 
until a minimum is reached for huj « 26 meV (20 meV) 
when the magnetic field intensity is 7 T (3 T). 

Comparing the top and bottom panels in Fig. |12[ it 
can be seen that the minimum (maximum) of the Fara- 
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day rotation angle coincides roughly with the maximum 
(minimum) of u'^^y. The latter fact agrees well with what 
could be concluded from the approximated result stated 
in Eq. ( 89 1 . In order to interpret the variation of the 



Faraday rotation angle with the photon energy, it is suf- 
ficient to use the simplified results derived in Sec. [IIP [ for 
T = 0, namely. Eqs. (57l-(63l. (This is clearly justified 
given the low temperature in the experiment in Ref. ,33. 
the respective thermal energy corresponds to about 0.01 
times the level spacing Ai = Ei—Eq [see Eq. (45 1 for the 
definition of A„] for both intensities of magnetic field.) 

For a magnetic field of 7 T (3 T), intraband transi- 
tions n = 9— ?'n = 10(n = 22— 7'n = 23) control 
the variation oi Op, from positive up to negative values, 
as the photon energy varies. Here, the index n denotes 
LLs with energy given by En = sigxilri) 2\n\hv p / 1 b [see 
Eq. ([5]) and text thereafter]. The remaining transitions 
contributing to the Hall conductivity are interband-like 
and occur at much higher photon energies huj ~ 'iEp, 
and thus it does not infiuence the Faraday rotation in 



the range of energy plotted in Fig. 12 

In this example, intraband transitions involve a very 
small difference in energy, even when the magnetic field 
is 7 T. The value of the intraband gap [Eq. (45 1] is Ajv^^ ~ 
16 meV (Aat^ ~ 7 meV) for B = 7 T (B = 3 T), which is 
comparable to F (here Np denotes the last occupied LL 
for a given Fermi energy). The exact calculation shows 
that the extrema points of the real part of the intraband 
Hall conductivity [Eq. ([ST])] occur at cj = 0, and. 



n 



A2 



r2 ± 2r 



r2. 



(92) 



Substituting the values given in the caption to Fig. [12] 
into the latter formula, we obtain lo™^''^ ~ 27 meV 
(^mtra ^ 20 mcV) for a field intensity of 7 T (3 T). As 
mentioned above, these are the points where the Faraday 
rotation reaches its minimum value. Increasing further 
the photon energy, fiiu > /iw"*''*, the Faraday rotation in- 
creases toward zero, essentially because at large u, below 
the interband threshold, the Hall conductivity becomes 
very low (Fig. 12 1 and no distinction arises between (t_ 
and (T+, and thus w t_. Increasing the photon energy 
up to hu! ^ 2Ep, the interband transition comes into play 
and drives the Faraday rotation. Interband transitions 
are important in samples with low electronic densities, 
as explained in the following section. 

The curves for 6p, computed either from the semiclas- 
sical expressions for the conductivity [Eqs. ( [73| a nd (74l] 
or via the EOM expressions [Eqs. (42 1 and (|55p] are al- 
most indistinguishable (see Fig. [12] top), in the range 
of photon energies considered, except for huj w 10 meV, 
where a very small deviation is observed when the inten- 
sity of the magnetic field is 7 T. 

The agreement between the quantum and the semi- 
classical solutions is explained by the similarity of the 
intraband gap A^r^^ and the cyclotron energy hujc [see 
Eq. ([tTJ]. The values for these quantities are AjVj^ — 
6.62(15.6) meV and huc ^ 6.58(15.4) meV for a field 
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FIG. 13: Low electronic density limit. Top: Faraday rotation 
angle (given in degrees) for free-standing graphene (e^ = 1) 
for different LL occupations: from left to right, Np ~ 0, 1, and 
2. The magnetic field intensity is B = 7 T, F = 10.5 meV and 
T — 0. Adding a dielectric substrate to graphene decreases 
the maximum amount of Faraday rotation that is achie vabl e, 
without introducing major qualitative changes [see Eq. (89l]. 
Bottom: Real part of the quantum conductivity tensor for 
the Fermi energies considered in the top panel. 



of 3(7) T. The agreement between the methods breaks 
down near the interband threshold, hu ~ 2Ep ~ 0.6 eV, 
where the quantum contribution arising from the inter- 
band transition cannot be neglected. 



C. Quantum jumps in the Faraday rotation: the 
low electronic density limit 

When low Fermi energies are considered, energy quan- 
tization becomes important (see Sees. II D and II E[ ). 
The limiting case occurs for < Ep 



< 



\/2hvp/li 



i.e., Np — 0. In this case, at T = 0, LLs with n > 1 
are all empty, and a single type of transition contributes 
to the Hall conductivity, n = — > n = 1. Since this tran- 
sition is interband-like, it cannot be explained within the 



semiclassical treatment (Sees. ]IID|and II El. This situa- 



tion is illustrated in Fig. 13 (bottom): when Np — 0, the 
real part of the Hall conductivity ]solid (green) line] has 
a finite (nonzero) value around u ~ {Ei + Eo)/fi. [Note 
that the extrema of the interband Hall conductivity can 
be obtained from Eq. ( 92 1 by making the replacement 
Anj, ^ hAflNf, with AflNp given by Eq. (|48l.] The 



Faraday rotation given by the semiclassical model is ob- 



viously ]dashed (red) line] since Ep ^0 ]Eq. (74l]. The 
respective Faraday rotation angle (top) is approximately 
proportional to —axy{u). 

At higher Fermi energies (i.e., Np > 0), two types of 
transitions contribute to the Hall conductivity: in gen- 
eral, for En„ < Ep < Emp+1, with Np > 1, the allowed 
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transitions are (i) interband between the hole's LLs with 
n = —Np and the electron's LLs with n = Np + 1 and 
(ii) intraband between LLs with n — Np and n = Np + 1 



(Sec. II D I. The maximum intensity of a' falls off with 



the inverse of the energy difference associated with a 
given electronic transition [Eqs. (62) and (63l]. Since, 
up to a good degree of approximation, the Faraday effect 
is controlled by a'^y, the latter means that the amount 
of Faraday rotation induced by the interband transitions 
at w = Ar^ATy will be smaller than the Faraday rotation 
due to intraband processes. 

The above-mentioned facts can be appreciated in 
where numerical data for 9p (top). 



Fig. 
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cr^j, and 

(bottom) are shown with Fermi energy increasing 
from left to right. As higher LLs in the conduction band 
become occupied, the spectral weight for the interband 
contribution to a'^y shifts toward higher energies (that 
is, hAfljyp increases — uj^}^^'^ increases). The opposite 
occurs for intraband transitions, since in this case, the 
relevant energy scale A^^ decreases with increasing Ep. 
As a result, the intraband part of a'^y concentrates its 
spectral weight at the lower edge of the plotted spectrum, 
and displays a much higher amplitude than its interband 
counterpart, as explained above. Similar conclusions ap- 
ply to 9p, as direct inspection of the bottom and top 
panels shows. 

When Np = 1, a significant departure from the 
semiclassical behavior can be appreciated in the intra- 
band region (0 — 100 meV). Remarkably, though, al- 
ready for Np — 2, the semiclassical Hall conductivity 
approximates the quantum result well, with a significant 
deviation occurring only near the interband threshold 
(«300 meV), where the semiclassical approach must nec- 
essarily fail. These features are in accordance with the 
general conclusions drawn in Sec. |IID[ 

For comparison, the real part of the longitudinal con- 
ductivity is also shown in the bottom panel in Fig. [13] 
The longitudinal current can be induced by photons 
which are resonant with any interband transition allowed 
by the Pauli principle (i.e., huj > hAflN^ ~ 2Ep), and 
hence many absorption peaks can be observed. On the 
contrary, the Faraday rotation essentially depends on a'^y 
and therefore is driven only by two resonances. 

Dependence on the Fermi energy and magnetic field — 
We have seen that when few LLs are occupied, quantum 
effects come into play and the semiclassical solution no 
longer gives an accurate description of the Faraday effect. 
The latter can even happen in the intraband region (see, 
e.g.. Fig. 13 middle panel), embodying the departure of 
the intraband gap A jv^ from its semiclassical analog, the 
cyclotron energy, hujc (see Table |l]). Given the importance 
of energy quantization for a low electronic density, we 
expect 0p to display abrupt behavior when the Fermi 
energy crosses the first few LLs. The latter behavior 
should reflect directly the step structure of the optical 
(or ac-) Hall conductivity a'^y{uj)^ 

Figure [M] shows the Faraday rotation angle versus Ep 
for a fixed magnetic field, B — 7 T. The heights of the 
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FIG. 14: Quantization of the Faraday effect in graphene. 
Top: Faraday rotation angle (in degrees) for free-standing 
graphene as a function of the Fermi energy at a magnetic 
field of B = 7 T for fio; = 10 meV (left) and huJ = 50 meV 
(right). The respective semi-classical result is plotted by the 
dashed lines. Other parameters: F = 10.5 meV and T = 12 K. 
Bottom: Same as the top panel but with F = 2 meV. 



steps are not uniform since the optical Hall conductiv- 
ity no longer obeys the dc quantization rule [Eq. (64l] 



When the Fermi energy crosses higher LLs, the smooth 
semiclassical result (dashed curves) is recovered. 



Combining the approximated formula for 6p [Eq. (89 1 

< 



valid for cfia'^luj) <C 2 and 6p 
conductivity at T = [Eqs. 



(57 



1] and the exact Hall 



and (58 1], explicit for- 
3e obtained. When the 
> 1, the expression 
Nevertheless, 



mulas for the step heights can 
Fermi energy crosses LLs with n 
for Adp becomes somewhat cumbersome 
simple analytical expressions can be obtained in some 
regimes. For instance, when the photon energy is small 
compared to relevant scales, /iw <C F <C the steps are 
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predicted to be approximately uniform, 



(n-i-n+l) 



4a ~ 0.03 rad, 



(93) 



where a denotes the fine structure constant, a ~ 
/ {ATrheoc) . In Ref. [19], for estimation of the magni- 
tude of the effect it was assumed that the step height of 
a'^y is approximately given by A<T'^y{uj) ~ /h, resulting 
in A9p ~ a. Rigorously, the step height for the transi- 
tions n = 0— >n = lis about 4e^//i, hence explaining 
the extra factor of 4 in our expression. In fact, in the 
limit fiw <C r <C -El, the steps in the Hall conductivity 
will all have approximately the same height, as in the dc 



case [see Eq. (64l] 



In Fig. [14] a decrease in the step's height relative to the 
estimated value in Eq. ( 93 ) can be observed already for 



the first step. This happens because the condition F <C 
El is too restrictive, and hence we relax this condition 
to F < but, at the same time, keep the low-photon- 
energy condition, hui <^ T. Doing so, leads to a better 
approximation , 



1 



1 + (6 + 471 + dnfl)j^ 1 + 4ri72 



X Aa , 



(94) 

where we have defined the dimensionless parameter 7 — 
T/Ei. Using this parameter, the validity condition of 
Eq. (|94| reads, 7 < 1 and fio; < F. 

Two physical scenarios where the Faraday steps are 
not imiform are shown in Fig. [T4| In the bottom panel, 
the transitions 71 = — > n = 1 {Ep ~ 100 meV) come 
with a variation of Op of roughly 1.8° (~ 0.031 rad) 
for hu! = 10 meV, versus —5.1° (~ —0.089 rad) for 
huj —50 meV, which does not agree either with the rough 



uniform estimative or with Eq. ( 94 1 . The reason for this 



discrepancy is that the condition hu <^T is not fulfilled 
for the photon frequencies considered in Fig. [14] Recall 
that in graphene, F is about about 10 meV, and thus 
infrared photons have fiw > F. It is therefore useful to 
derive approximate formulas for A9p that are valid in 
the regime huj ^T. Defining uj — Ei/{huj), we arrive at 



A9 



>n+l) 



1 



1 - 2(1 + 2n)w2 + ^4 1 _ 2(3 + 2n)J)2 + Cb^ 

(95) 

Substituting for the respective values of w, we obtain 

Ae'^p^^'' = 1.8° and A6'^"^^' = -5.3°, for huj = 10 eV 
and huj — 50 meV, respectively, which agrees well with 
the numerical results reported in Fig. [14] for F = 2 meV. 
As for the steps observed in the top panel in Fig. [14] 



they cannot be explained accurately with Eq. ( 95 1 since 



in that case we have huj « 0{T). We stress that Eqs. (94) 
and (95) are only accurate when the statement, Eq. (89), 



provides a good description of the Faraday effect in 
graphene, which in practice means very high photon en- 
ergies hu) (see also Fig. 111. For the parameters used in 



analytical expressions for A9p are accurate only for the 
first few steps. 

Figure [Ts] shows the variation oi 9p with the magnetic 
field for two cases; (i) low doping (Ep — 0.05 eV) and (ii) 
high doping {Ep = 0.3 eV). In the latter case, we are well 
inside the semiclassical regime even for the maximum in- 
tensity of the magnetic field considered {B =7 T), and 
thus no distinction can be made between the curves com- 
puted using the semiclassical conductivity tensor or the 
EOM formulas. In this regime, the Faraday effect in- 
creases monotonously with the magnetic field. 

For a low electronic density, on the other hand, the 
agreement between the Boltzmann and the EOM for- 
malisms takes place only for low magnetic fields. For 
increasing values of the magnetic fields, such agreement 
ceases to occur as soon as the intraband gap does not 
match the cyclotron energy huJc- Then, energy level 
quantization becomes important and the EOM expres- 
sions must be considered (i.e., Np is small; see Sec. HE): 



this explains the departure from the semiclassical value 
for 9p observed in the right panel in B x 1 T for 
hoj = 10 meV {B w 0.5 T for = 30 meV). If the 
magnetic field intensity is higher than a given value, we 
necessarily have Np = (for 0.05 eV this value is about 
1.9 T). In this case, the Hall conductivity, at T = 0, is 
fully determined by a single type of interband transition, 
and, assuming Ei[B) ^ /iw,F, we obtain [see Eq. (58)], 



2e2 
h 



2a ~ 3 X 10" 



(96) 



The latter considerations explain the plateau formed at 
_B w 2 T [solid (blue) line] in the right panel in Fig. [15] 
The dashed-double-dotted (red) line corresponds to pho- 
tons with a higher energy, shifting the formation of the 
plateau toward higher fields. Equation (96) is indeed the 



Fig. [14| where the photon energies are not very high, our 



high-magnetic-field limit \Ei{B) ^ energy scales] of the 
Faraday rotation induced by single-layer graphene. 

Although the measured Faraday rotation angle is re- 
markably large given that it comes from a single graphene 
layer, in both low and high doping regimes (see Fig. 15), 
it is the needed magnetic field, i? > 1 T. The goal is then 
to obtain large Faraday rotation angles using graphene 
and modest fields at the same time. A simple idea that 
uses the nonreciprocity of the Faraday effect is to enclose 
graphene between two mirrors. We discuss this possibil- 
ity in the following section. 



D. Enhancement of Faraday rotation in a cavity 
geometry 

We have seen that the existence of intraband and inter- 
band transitions in graphene permits the generation of fi- 
nite (nonzero) Faraday rotations in different ranges of the 
electromagnetic spectrum. In doped graphene {Np > 1), 
e.g., the intraband gap is bounded from above by 



E" 



<Ai^E2-Ei~ 15/BmeV • T"^/2 , (97) 
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FIG. 15: Faraday rotation angle (in degrees) as a function of 
the magnetic field for Ef = 0.30 eV (left) and Ef = 0.05 eV 
(right). In each panel two photon energies are represented — 
hijj — 10 meV [solid (blue) line] and fiw = 30 meV [dashed- 
double-dotted (red) line] — with the respective semiclassical 
counterparts shown by dashed lines. Other parameters as in 
the top panel in Fig. |14| 



implying that, by using magnetic field intensities ^ 1 T, 
graphene can be exploited for magneto-optical applica- 
tions from the microwave up to the far-infrared regimes 
/ = E/h < 3.6 THz (an example of terahertz Faraday ro- 
tation driven by intraband transitions is given in Fig. 12 I. 



Another possibility is to make use of transitions connect- 
ing the valence and conduction Dirac cones, whose inter- 
band gaps are bounded from below. 



E" 



> hA^i ^Eic^ 36 VB meV • T 



-1/2 



(98) 



thus obtaining far-infrared up to visible light light fre- 
quencies (an example of mid-infrared Faraday rotation 
driven by interband transitions is shown in Fig. 13 1. We 



recall that increasing the electronic density in order to ob- 
tain even larger interband gaps (/iAr2„ with n > 1), and 
thus shifting the magneto-optical response of graphene 
above the midinfrared, huj ^ hAflNp — '^Ep, creates op- 
tical Hall conductivity peaks with a low intensity. As a 
consequence, very small Faraday rotations are produced 
already in the near- infrared regime. A good estimate for 
the maximum achievable interband-induced Faraday ro- 



tation can be obtained from Eqs. (63 1 and (89 1 



max \9p\ 



V 2loV 



X a , 



(99) 



which, for example, taking i? = 7 T , F = 10 meV, 
and hw = 1 eV, leads to max 1 0^1 — lO^'^. Although 
the amounts of terahertz Faraday rotation, hw ~ OmeV, 
reported in our figures are well within state-of-the-art 
capabilities [the resolution for Faraday measurements in 
terahertz time-domain spectroscopy is presently limited 



to 1 mra(P^('^0.06 degrees)], high magnetic fields, ~ 1 T, 
are still needed which can be a disadvantage for specific 



applications; moreover, according to Eq. ( 99 1 the needed 



magnetic field increases as higher photon frequencies are 
to be probed. 

The situation is very different in other 2D electron 
gases, for which Op is proportional to the sample's thick- 
ness (as the light travels farther through the mate- 
rial, more Faraday rotation accumulates). Single-layer 
graphene, on the other hand, being one atom thick and 
hence truly 2D, requires the use of high magnetic fields in 
order to detect Faraday rotations. It is therefore natural 
to ask whether it is possible to conceive a setup leading to 
accumulation of Faraday effect; ideally, such setup would 
avoid the use of several samples and, at the same time, 
take advantage of the broad magneto-optical response of 
single-layer graphene. 

In what follows, we discuss a graphene-based system 
that can enhance the intrinsic graphene's Faraday rota- 
tion at any frequency and thus can cope with the diffi- 
culty mentioned above. The idea consists in enclosing 
graphene in an optical cavity: due to intracavity in- 
terference, photons undergo several round trips within 
the cavity before leaking out. Loosely speaking, due 
to nonreciprocity of the Faraday effect, accumulation of 
Op then takes place each time a photon passes through 
graphene; a sketch of the experimental apparatus is 
shown in Fig. [TH) 

Explicit calculations (see below) show that giant Fara- 
day rotations are achieved even when the optical finesse 
of the cavity is modest. The optical finesse can be eas- 
ily tuned by changing the reflectivity of the end mirrors: 
the higher the latter quantity, the larger is the number 
of round trips of photons inside the cavity, and hence 
further Faraday accumulation occurs. Indeed, the cav- 
ity geometry gives a straightforward solution to mimic 
the effect of a sample's thickness (absent in single-layer 
graphene) . 

Following the steps in Sec. [Ill A] we write the boundary 
conditions of the electromagnetic field in terms of circu- 
larly polarized waves. Employing similar notation, we 
define the input and output circular vector amplitudes. 



5T* 



{E^,El) 



(100) 
(101) 

(the reflected and transmitted waves £'j_ having anal- 
ogous deflnitions). The flrst (second) component of the 
vectors, Eq. (100 I and (101 1, refers to the complex am- 



respectively (see also Fig. 16 1 



where E'l = E"^ 



plitude of light traveling in the positive (negative) z di- 
rection. 

The output field, 5™', and thus the total Faraday ro- 
tation angle, can be more conveniently computed using 
the transfer matrix formalism. The method is explained 
in detail in Appendix Here, we just state the basic 
results: the T matrix, by definition, connects the input 
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output field 



L 



FIG. 16: Schematic of the graphene-optical cavity system: 
linearly polarized hght shines into an optical cavity with 
graphene placed at the center. The field inside the cavity 
perceives graphene as an extra boundary and hence the two 
halves of the cavity operate as independent cavities of effective 
size L/2. Matching the light frequency hui with a resonant fre- 
quency of the cavity cj = rnrc/L (n £ N ) traps photons inside 
the cavity for several round trips. As a consequence, Fara- 
day rotation accumulates due to multiple passages through 
graphene, leading to an output field with a large Faraday ro- 
tation. 



and output vector amplitudes, according to 



(102) 



where is a product of individual transfer matrices 

for each boundary (optical component, metallic surface, 
etc.). Its inverse permits us to compute given the 

input field S^^, and hence the optical characteristics of 
the cavity-graphene system. In particular, the circular 
transmitted amplitudes, t± = E^./. 

For the geometry posed in Fig. 
T matrix reads 



[J]^ , are given by t± = 
16] the input-output 



T" 



in 



jL/2c q 







■Tl 



-iujhjic Q 







^iujL /2c 



Each operator in Eq. ( 103 1 propagates the electric field 
to the right until a boundary is reached. encodes the 
effect of the first interface, a mirror, and depends only 
on the mirror's transmission and reflection amplitudes, t 
and r, respectively. It can be written as 



1 



1 



'I 



(104) 



(For a derivation, see, e.g., Ref. [SIl) After interac- 
tion with the left-end mirror, photons can enter into the 
cavity and propagate for a distance of L/2 before the 
next interaction. This means that another T matrix is 
needed; free propagation merely adds a phase to the dec- 
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FIG. 17: Faraday rotation angle of a cavity-graphene system 
in the semiclassical and quantum regimes. Top left:^^ as a 
function of the photon energy for a cavity-graphene system in 
a magnetic field of 7 T . The Fermi energy reads Ep = 0.3 eV 
and the cavity mirrors have r = 0.99. Other parameters: 
r = 10.5 meV and T = 12 K. Top right: 0f versus the 
refiection amplitude r for Hlj = 19 meV. Bottom: In the left 
Tm ■ (103) (right) panel, the Fermi energy reads Ef = 0.05 eV {Ef = 
0.1 eV), which corresponds to an LL occupation of Nf = 
(Nf = 1). The orange dashed line shows 6i? as obtained with 
the semiclassical conductivity tensor. 



represented by a diagonal matrix, which is the second 
operator in Eq. (103 1. At z = L/2, photons arrive at 



trie field [see Eq. (A5l and text thereafter] and thus is 



the air-graphene-air interface, whose T matrix we de- 
note by T^. (More involved types of interfaces could 
be considered: for example, air-substrate-graphene-air. 
The present choice has the advantage of keeping the 
mathematical expressions elegant; generalization to other 
configurations using the present formalism is straightfor- 
ward.) The graphene's T matrix depends on the mag- 
netic field intensity, electronic density, temperature, and 
LL broadening, via the complex optical conductivity of 
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graphene a± (lu) ; its explicit form is 



Zo(T=p(w) 

2 - Zo(T=p(a;) 



(105) 



where Zq — iiqc denotes the vacuum impedance; see Ap- 
pendix 1^ for a detailed derivation. Finally, the second 
line of Eq. ( |103[ ) propagates the field in free space for a 
distance of L/2 and adds the right-end mirror. 

The Faraday rotation angle is obtained from 9f = 
(l/2)arg(t+/i_), with the circular amplitude ratio t+/t^ 
given by [Ti"^°"*]i,i/[r|"^™']i,i [Eq. (IAT])]. After some 
algebra, we arrive at 



t+ ^ 2 + Zpa+juj) - |r| [Zpa+juj) - 2] e""^/^ 
t_ 2 + Zoo-_(w) - |r| [Zo(T_(a;) - 2] 6*"'^^/'= 



(106) 



from which 9p can be immediately deduced. Setting r — 
in the latter expression leads to the previous result in 
the absence of a cavity [compare with t± as obtained from 
Eq. (|83]) with er = l]. 

When r > 0, interference takes place and photons 
can make several round trips before being transmitted 
through the cavity. On an intuitive basis, we then ex- 
pect that the Faraday rotation angle can be enhanced 
due to multiple passages of photons through graphene, 
which indeed is the case, as shown in Fig. [TT) Hereafter, 
the size of the cavity is set to L = mrc/uj, with n odd. 
The solid line shows Op for the cavity-graphene system 
and the dashed-dotted line shows 0f for free-standing 
graphene for the same parameters: clearly, in the range 
of frequencies considered, the Faraday effect is greatly en- 
hanced. For example, for a low frequency, hio w 10 meV, 
9f has increased by a factor of about 5, reaching a value 
of 55 degrees, whereas for Huj ~ 19 meV, Op increases by 
a factor of about 20, reaching a value of approximately 
25 degrees. 

Direct inspection of Eq. (106) discloses the observed 
boost of the Faraday effect: when r — > 1 and the phase 
factor exp(ja;L/c) = — 1, the constant factor of 2 cancels 
in both the denominator and the numerator, leading to, 

, cr+(cj) 



2ief 



n odd (7_ {uj) ' 



(107) 



which can present large arguments, 26 p. The opposite 
limit, r — 0, in which the isolated graphene system is 
recovered, leads to much smaller arguments, since gen- 
erally 2 3> Z{)lTaij±, which implies that the real part of 
Eq. (83 1 is predominant. Choosing a cavity mode with n 
odd and r ~ 1 is fully equivalent to taking a large num- 
ber of equally prepared graphene sheets placed in a row 
(Appendix [bJ . The cavity geometry therefore permits us 
to take advantage of large Faraday rotation accumulation 
using a single graphene sheet. 

In a cavity geometry, the Faraday rotation is no longer 
dominated by the behavior of (t^j^(w) [see Eq. (89)[, for Op 
now depends on the full conductivity tensor [Eq. (107l[. 
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FIG. 18: Faraday rotation boost in the infrared and visible 
ranges. Left: Faraday rotation angle versus photon energy 
of a cavity-graphene system with Ep ~ 0.85 eV (top) and 
Ef = 0.3 eV (bottom). Right: Transmissivity of a cavity- 
graphene system for the same parameters considered at the 
left panels. Inset: Transmissivity of intrinsic graphene for the 
same parameters. Other parameters as in Fig. |17[ In order to 
obtain a measurable Faraday rotation at ticj « 1.7 eV [solid 
(red) line[ [huj « 0.5 eV (infrared)[, it is necessary to tune the 
intraband resonance according to Ep — hui/2. 



The most visible consequence of the latter fact is that 
photons with htj w 20 meV undergo considerable Faraday 
rotation angles in a cavity geometry, whereas, in a sin- 
gle passage through graphene, photons with such energy 
do not produce Faraday rotation at all (Fig. 17 1. This 
apparently counter-intuitive result is due to induced ellip- 
ticity in single passages and is explained in Appendix [B] 

a. Semi-classical versus quantum regimes in a cavity 
geometry — Figure [Tt] (top) considers the case of Ep = 
0.3 eV and B = 7 T, well inside the semiclassical regime, 
for which the ac conductivity is dominated by intraband 
contributions over a wide range of frequencies (Sec. [ll| ; 
the corresponding intraband Faraday rotation is seen to 
be greatly enhanced in the cavity geometry. 

The low-electronic-density regime of the cavity- 
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graphene system is shown in the bottom panel in Fig. [17] 
Remarkably, for energies above the interband thresh- 
old, namely, hw > Ei ~ 95 meV for Np ^ (left) 
and hu > El + E2 230 meV for Np ^ 1 (right), 
9p{ui) presents a behavior qualitative different from that 
of an isolated graphene sheet (black dot double-dashed 
curve): oscillations do emerge. These oscillations are hin- 
dered in single-photon passages through graphene (see 
also Fig. 13 1, but for multiple-photon passages, in the 



high-frequency limit, Shubnikov-de Haas oscillations in 
the longitudinal conductivity axx{'^) (Fig. |4| are criti- 
cal in defining the orientation of light polarization axes. 
These oscillations are obviously absent in the semiclassi- 
cal Boltzmann calculation [dashed (orange) curve]. In the 
top panel, where Ep = 0.3 eV, such oscillations are not 
present because the represented photon energies are well 
below the threshold for interband transitions huj ~ 2Ep. 

h. Near- infrared and visible-range Faraday rotation — 
We finish this section by mentioning an important ap- 
plication of the cavity-graphene system: interband- 
induced Faraday rotations in the near-infrared and visi- 
ble regimes. Figure [T8[ shows that energetic photons can 
attain Op ^ 1 hy tuning the Fermi energy to sufficiently 
high values. In this regard, the top panel shows numer- 
ical data for graphene with Ep = 0.85 eV; such a high 
doping level of graphene samples is feasible using chem- 
ically synthesized graphene with ferroelectric substrates 
(instead of the conventional Si02).'^ 

Given the mirror reflection amplitude considered, r = 
0.99, photons are trapped for a large number of round 
trips. This means that it is highly probable that photons 
get absorbed by graphene before leaking into the cav- 
ity. This explains why the transmissivity of the cavity- 
graphene system, as shown in the right panel in Fig. [18] 
is well below 1 (but still large enough that the effect can 
be measured). One way of increasing the transmissivity 
of a cavity-graphene system is to decrease the quality of 
the mirrors, at the expense of decreasing the maximum 
achievable 9p. 

We, finally, remark that the nonlinearity associated 
with next-neighbor hopping t' in a honeycomb graphene 
lattice can play a role for photons with huj > 1 eV, 
and hence corrections to the Dirac cone approximation 
(Sec. II A), and thus to the EOM solutions, may exist; 
such corrections are expected to be very small, however.-* 



IV. CONCLUSION AND OUTLOOK 

In the first part of this work, the EOM method has 
been adapted to the study of magneto-optical trans- 
port of electronic systems. To illustrate the method, 
the magneto-optical conductivity tensor of single-layer 
graphene in the Dirac cone approximation has been 
derived, accounting for both intraband (semi-classical) 
transitions and interband transitions between the valence 
and the conduction bands. 

The general regularization procedure to obtain the reg- 



ular conductivity tensor from the solutions of the EOM 
for the current operator has been established; this pro- 
cedure is shown to lead to the correct formulas with- 
out the need for evaluation of the Kubo formula. To 
the best of the authors' knowledge, such a procedure has 
not been discussed in the literature so far. In addition, 
quantitative comparisons between the quantum EOM so- 
lutions and the semiclassical Boltzmann formulas, in the 
full optical spectrum, and in both low- and high-doped 
graphene samples, have been given throughout. 

In the second part, the Faraday rotation effect in 
single-layer graphene has been studied in detail; in partic- 
ular, simple formulas for the step heights in the quantum 
Hall regime have been derived. Our results have been 
shown to account well for available experimental data in 
the semiclassical regime. 

Finally, we have proposed a simple experimental appa- 
ratus based on an optical cavity that leads to an enhance- 
ment of the Faraday rotation effect of graphene by orders 
of magnitude, thus allowing to obtain giant Faraday ro- 
tation angles in the infrared region and modest Faraday 
rotation angles in the visible region. 

We hope that the present work further stimulates the 
research on magneto-optical properties of ultrathin 2D 
gases and graphene-based soHd-state devices. 
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Appendix A: Transfer matrix formalism 

The transfer matrix (T matrix) approach is a widely 
used method in optics and related fields and provides an 
efficient means of calculating the amplitude and phase of 
transmitted electric fields through an arbitrary number 
of interfaces. In this appendix, we give a self-contained 
review of the method and derive explicitly the T matrix 
for a general 2D conducting media. 



1. General formalism 

For concreteness, we assume that an incident electro- 
magnetic wave of frequency oj, travels in the z direction 
through a set of N metallic interfaces, placed normal 
to the direction of propagation, with labels a„_„_|_i, and 
located at positions 2; = z„ (n G 1,2,. ..TV). These in- 
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FIG. 19: Schematic of an optical system consisting of an array 
of interfaces separated by different types of dielectric media. 
An electromagnetic wave, coming from a medium 

with dielectric permittivity ei interacts with an interface Q12. 
As a result, it is partially reflected and partially transmit- 
ted into the medium £2. Equivalent events take place at the 
remaining interfaces. The vectors with superscript +(— ) de- 
note the component of the electric field traveling in the posi- 
tive (negative) direction of z. A uniform static magnetic fleld 
B = _Bey is assumed. 



terfaces are separated by dielectric mediums — Fig. 19 
shows the configuration we have in mind. 

The electric field is separated according to the direction 
of propagation: represents the part the electric 

field traveling in the positive direction of z, within the 
region n, whereas E~ [z) represents the part traveling in 
the opposite direction. 

As shown below, the calculation of transmitted and re- 
fiected amplitudes becomes easier by writing the bound- 
ary conditions in terms of circularly polarized waves (see 
also Sec 



plitudes. 



III A). Therefore, we focus on the circular am- 



EtA^) = Et,Az) + iTEt^^y, 



(Al) 



where t is the polarization index: r = ±1 [-|-1(— 1) means 
right-handed (left-handed) circular polarization] . Indeed, 
in a given region i, the total (complex) electric field is the 
sum of both components. 



E,,{z,t) = E+{z)e 



(A2) 



The physical electric field is obtained by taking the 
real part of the latter expression. We omit the time- 
dependence in the remainder of the appendix. 

The T matrix connects the amplitude of the electric 
field to the left and to the right of a given boundary 
(interface). Take for instance, the interface labeled ai 2 
in Fig. 19 The respective T matrix, T^'^, is defined as. 



E2M) 



1,2 ( E+Jz^) 
El Ml) 



(A3) 



where z^ denote the position where the electric field is to 
be evaluated: right after (+) or before (— ) the interface 
located at z = zi. For ease of notation, we define the 
vector of amplitudes. 



, r ( ^ ) 



E+A^) 

E-M) 



(A4) 



and drop the superscripts in the coordinates z„. 

If more than one interface is present, the light propa- 
gates a given distance before interacting with the next 
component. Propagation of light through a dielectric 
medium merely adds a phase to each r circular com- 
ponent of the electric field. Indeed, its action can be 
represented by a diagonal matrix. 









^iknAZn 



£n,r{zi+l) , (A5) 



where the index n just takes the values for which there is 
intermediate light propagation, i.e., n = 1, A'^ — 1, the 
wave vector depends on the dielectric medium according 
to fc„ = ujy/en/c and Az„ = z„+i — z„ is the width of the 
region n. Note that Eq. (A5l defines a particular case of 



a T matrix, which we denote by Fn- 

The problem of finding how the output electric field, 
of definite polarization r, immediately after leaving the 
last interface, E'j^j^^ t^zn), relates to the incoming elec- 
tric field, with the same polarization t, E^^{zi), then 
amounts to take the product of the individual transfer 
matrices, 



£i.t{zi) 



rf,l,2 £n rf,2.3 rfiN-l.N-fp rj^N.N+lp \ 

J^T ^2^r ■■■^T i^N-Lr +i,r{Zi^ ) . 



(A6) 

The total T matrix J'ln-^out ^^^^ ^^1^ desired information: 



E 



N+1,T 



Et 



1/ 



T.. 



in— >-out 



1,1 



(A7) 



As for the relation between the output field and the re- 
flected fleld at the flrst boundary, E^^i), we obtain. 



E 



N+1,T 



= 1/ 



^11 



(A8) 



In what follows, we show how to construct the T ma- 
trix for a general 2D conducting medium. Knowledge 
of the transfer matrices allows us to determine the char- 
acteristics of transmitted and reflected light through a 
general set of conducting 2D thin fllms, mirrors, etc., by 
employing Eq. (A6I. 



2. T matrix for a general conducting 2D interface 

We restrict the present derivation to non-magnetic me- 
dia, and assume the standard constitutive relations hold. 



J(r,cj) = (T(r,w)£;(r,a;) 



(A9) 
(AlO) 



where D, e, and a denote the displacement fleld, permit- 
tivity and conductivity, respectively. Also, and without 
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prejudice, we take the 2D conducting interface to be lo- 
cated at z — 0. The T matrix is defined as 



(afc) 



E+^{z = 0^ 



(All) 



where a (b) is the bookkeeping index for the medium at 
the left (right) of the interface. 

Various constraints emerge due to continuity of E (and 
its derivative) at the 2D conducting interface. Indeed, 
Maxwell equations imply that. 



Ea{0) = E,iO), 



(A12) 



fdEa 



dEt, 

dz 



where the conductivity tensor reads 
a = aij{uj)6{z) . 



iwMo<T-Eh(0) , (A13) 



(A14) 



The conductivity depends on the light frequency w, and, 
generally, also on other quantities (Fermi energy of the 
interface, temperature, etc.). In the latter expression, 
the subscripts i,j = x,y are Cartesian coordinates. In 
terms of circularly polarized fields, Eq. (A13l reads. 



ka{El, - E-,) - h{E+^ - E^^^) = c^AtoX 

X {E+, + E-^)<j_A^), 
(A15) 

where we have admitted an isotropic medium, axx = <^yy, 
and have defined 



(A16) 



The statement, Eq. (A15l, shows that the two circularly 



polarizations are decoupled, even in the presence of a 
complex conductivity a±{uj). This is why it is advanta- 
geous to write the boundary conditions in terms of cir- 



cularly polarized fields (Sec. Ill A) 



According to the definition of the T matrix [Eq. ( All l] 
we need to relate with E^^ and E^^ with E^^, sep- 
arately. To do so, we make use of the continuity condi- 
tion, Eq. (A12l, written in circular waves, E'^^ + E^^ - 



E, 



Ef^ ^, in order to arrive at. 



±2kaE^^, = hiE+^-E-^) 



[iJ^10C7^riuj)±kaKE+^+E^^). (A17) 



Combining Eq. ( All ) and the latter expression, we ar- 



rive at the desired result, 

T 

where 



1 



\ab \ab 



b 

T, 



A"' 



b 

r, + - 



A afc 



ka±kb± a;/zocr_r(w) . 



(A18) 



(A19) 



3. Example: T matrix of suspended graphene 

The T matrix of suspended graphene can be obtained 
immediately from Eq. (A16l. Admitting that the medi- 



ums at the left and right of the single-layer graphene 
sheet are air, we obtain 



jigraph 



2 + Zoa^'^^'^iiu) 

7^;^P^(w) 2 - ZaCT^T'\^) 

(A20) 

where Zq = fJ^oc is the vacuum impedance. 



graph / 



Appendix B: Faraday Effect 

In the present appendix, we derive the exact analytical 
conditions for the existence of Faraday rotation and dis- 
cuss their modification when graphene is enclosed in an 
optical cavity. Despite the focus on graphene, most of the 
conclusions drawn here apply generally for systems pos- 
sessing in-plane symmetry. Once again, for simplicity, we 
consider the case of suspended graphene; generalization 
to the case of graphene on top of a substrate is straight- 
forward using the general formulas given in Appendix [X] 



1. Conditions for Faraday effect in free space 

We consider a target graphene sheet, placed on the xy 
plane, subjected to a normally incident electromagnetic 
wave, linearly polarized along the x axis, ExC'^^'^*' . The 
magneto-optical Faraday effect takes place when a mag- 
netic field B = Bcz is applied. Then, Lorentz force acts 
on free carriers, producing a Hall electronic ac current, 
which, under specific conditions (see below), will produce 
out-of-phase radiation polarized transversely to the im- 
pinging field, EyC^^'^'^e^'^ . As a consequence, the resulting 
electromagnetic wave sees its polarization plane rotated. 

Without loss of generality, consider the graphene sheet 
to be placed at 2 = 0. In the circular basis, Er = 
(l/2)(ea; -I- TiEy), the electromagnetic field aX z — 0~, 
reads 



E{Q-) = Eoe-'''^\e+ + e_) . 



(Bl) 



Note that the actual electric field is given by the 
real part of the latter equation. After interaction with 
graphene, each of the circular components t — ±1 
changes according to Eq. (A20). The field right after 



the graphene plane is given by 



E{Q^ 



Eoe- 



1 + /3- 



1 + 



(B2) 



with I3± = Zocr±(w)/2. To determine whether the plane 
of polarization has rotated, we write the latter equation 
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in the Cartesian basis, 



E{0+) 



Eoe- 



2(l + /3+)(l + /3_) 
X [(2 + Zaaxx) ex — Z^Uxyey] , 



(B3) 



where we have used the definition of I3± to sinipHfy the 
term inside brackets. Obviously, no Faraday rotation 
takes place when cTxy{uj) = 0. On the other hand, hav- 
ing axy{^jj) 7^ does not suflice to rotate the polarization 
plane; linear polarization can change to elliptic polariza- 
tion with main axes along x and y (this is the case for 
_B = 5 T and huj « 15 meV, as shown in the top panel in 
Fig.[TT] elliptic polarized light leaves the graphene sheet, 
5 ~ 0.15, but still Op = Q)- For this reason, the actual 
condition for the existence of Faraday rotation is 



\axy \ > A Arg 



° ^ ^±(2m + l)-,meNo, 



ZqCTx 



(B4) 

The amount of Faraday rotation is given by Eq. ( 85 ) and 



thus can be obtained directly from Eq. ( B2 1 , reading 



?F = ^Arg 



Zna_ 



Zoa+ 



(B5) 



In many situations (e.g. high photon energies and high 
electronic density), the longitudinal conductivity obeys 
Zoa'r^^ <C 2 -|- Zqct^j., thus leading to the approximate 
condition, \a'^y\ > ^ 9f > 0. This is consistent with 
the approximated formula derived for the Faraday rota- 
tion angle [Eq. (89l], which states that 9p is proportional 



to a' (see also Fig. Ill 



Conditions for the Faraday effect in an optical 
cavity 



In Sec. |IIID| we have seen that large Faraday rotations 
Op can be achieved in a cavity-graphene system, even for 
such photon energies that do not cause Faraday rotation 
in free space. An example is given in Fig. |17| in free 
space, impinging light with huj w 20 meV does not change 
its polarization direction, 6p ^ 0, whereas Op can be as 
large as 25° for graphene mounted on a cavity geometry. 

In order to explain the above-described phenomenon, 
it is sufficient to consider the simplified situation where a 
normally incident photon interacts with graphene twice 
in a row. For concreteness, we take two graphene sam- 
ples, equally prepared, separated by a given distance W. 
Let the photon frequency ui be such that no Faraday rota- 
tion is produced in the passage through the first graphene 
sample, that is. 



Arg 



2 + Zoaxxi^) 
Zo(Jxy{u)) 



±(2m + 1) 



2 ' 



(B6) 



for some m € Nq [see Eq. (B4|[. In the latter expression, 
it is assumed that axy{(jj) 7^ which is the case when a 



magnetic field is present. In these conditions, after the 



first passage, the electric field [Eq. (B3l[, can be written 



El 



Eoe- 



2[l + /3+(LD)][l + /3„(tS)] 
X [|2 + Zocrxx{i^)\ ex ± i \Zoaxy{oj)\ ey] 



(B7) 



where (f> — Arg[2 -I- ZoCFxxiijj)] and the sign ± depends on 
the actual argument of axy{ui). The latter equation de- 
scribes a field elliptically polarized with main axes along 
X and y (i.e.. Op = 0). We thus see that although no 
Faraday rotation occurs when Eq. (B6) is fulfilled, the 



polarization changes from linear to elliptic, an unavoid- 
able consequence for Lorentz force enforces some radia- 
tion to be emitted that is polarized along the y axis. 

In order to determine the field after the second passage, 
and hence demonstrate our point, i.e., that some Faraday 
rotation must necessarily be produced in multiple pas- 
sages through graphene (such as in a cavity geometry), 
we make use of the transfer matrix formalism. Indeed, 
we approximate the total T matrix by f graph . j|;raph ^^^lis 
approximation is exact when the phase for free propaga- 
tion between the graphene sheets, loW/c, equals 2TO7r). 
Employing Eq. (A7l, we obtain 



Eoe- 



[l + 2/3+(cS)][l + 2/3_(cS)] 
X {[1 + Zoaxx{i^)] {u)ey} . (B8) 



This time, the condition for zero Faraday rotation. 



Arg 



1 + ZoCTxxi^) 



Zoaxyi^}) 



= ±(2to + 1)^ 



(B9) 



cannot be fulfilled because Eq. ( B6 1 fixes the photon fre 



quency in this example. Then a finite (nonzero) Faraday 
rotation is produced in the second passage. 

The case of graphene in a cavity geometry is more in- 
volved because intracavity interference takes place. Nev- 
ertheless, the physics behind the boost of Faraday rota- 
tion is analogous: if, for graphene subjected to a trans- 
verse magnetic field, it turns out that the first photon 
passage yields Op = 0, then, in the following passages it 
must be that 9p > 0. See, for instance, Eq. (1071, valid 



for an optical cavity made of mirrors with very high re- 
fiection amplitudes: because |cra;y(a;)| > 0, for B > 0, 
then Op > for all light frequencies. 



3. Row of graphene sheets 

Taking a number N of graphene sheets separated by 
W, such that ujW/ c = 2imr, leads to the following electric 
field, right after the last graphene plane: 



En — Eqc 



1 



1 



1 + 



1 + N[3+ 



, (BIO) 



26 



and hence in the limit iV ^ 1 we obtain, 



g+(^) 
(T_(a;)' 



(Bll) 



which coincides with the result obtained for the cavity- 
graphene system, given by Eq. (107 1. 



Appendix C: Regularization of the EOM optical 
conductivity 

The EOM approach consists in extracting the opti- 
cal conductivity from the average of the current opera- 
tor J{t) (obtained through the corresponding Heisenberg 
equation) . 

This method avoids the calculation of current correla- 
tions, and hence short-circuits the calculation of (Tij(u)). 
The crucial point of the EOM approach is the regular- 
ization of the following expression. 



where 0{uj) {O = J, E) is defined via 
0{t) = 0(cj)e-*"* +C.C 



(CI) 



(C2) 



Equation. \C2\ is valid for a monochromatic electromag- 
netic field, A = Aoe"^* -I- c.c, and for EOM solutions 
Jiiuj) in first order in Aq. For convenience, we write 
the external electric field as E(t) = E+{t) + E^{t), with 
E±{t) = ±iujAoe'f''^^. 

Despite the resemblance of Eq. ( |C1[ ) to the Ohm's law, 
tjjij^ijj) is not the optical conductivity: in the linear re- 
sponse regime, the EOM solution can be put into the 
form 



Jit) = ilj{Lu)E+it) + c.c. 



(C3) 



with -0(0;) as defined in Eq. (CI I. On the other hand 



the conductivity, a^t), is defined via the relation 



Jit) 



dTa{t - t)E(t) . 



(C4) 



The Fourier transform of Eq. ( C4 1 is nothing more 
than Ohm's law, J{Ld) — a{uj)E{uj), with (j{uj) = 
J^^dte''''^+'°^^'^cr{t). The function a{uj) is analytic in 
the upper complex plane and therefore satisfies Kramers- 
Kronig causality relations. 

From Eqs. (C3l and (C4), we immediately conclude 
that, '0(w) ^ ct(cij). The bottom line of the EOM ap- 
proach is that the tensor tp{uj) can be exactly transformed 
into a'(ijj) via a simple regularization procedure, as we 
show in what follows. 

Without loss of generality let Aq = Aoe^, with Ao e M, 
and consider that no current fiows in the absence of ex- 
ternal perturbations, (J^jit)) =0. Since we are interested 



in the regular part of the optical response, we also take 
J^jit) = Jj' (t) = Jj{t); then, in first order in ^o, 



{Mt))H 



drA(r)([ji(T),j/(i)])0, (C5) 



with i — x,y. Using the Lehman representation, and 
similar notation as employed above, the latter expression 
can be written as 



{Mt))H 



zh 



X e 



E 

(r-t) 



c?rA(r)(m| Jj.|n) (n| Jj|m) x 



-I3E„ 



(C6) 



Since we wish to find the explicit form of ip{uj), we per- 
form the integration over the variable r. We obtain 



{Mt))H 



-y 



1 



Z/i ^ a; -f a;„m 



X (e-/^^" 



^0^ 



{m\Jx\n){n\Ji\m)x 



c.c. 



(C7) 



where a small imaginary part has been added to ensure 
convergence. Making use of the definition, Eq. (C3l, we 
arrive at the desired result. 



1p^J{uJ) 



1 1 

Zh iuj 



E 



1 



X {m\Jj\n){n\J,\m) {e'^^" - e 



-0E„ 



(C8) 



where i — x. We also have ^xy{^) = —"^Pyxi^)- 

On the other hand, the frequency-dependent conduc- 
tivity is obtained from the Fourier transform of cr(i), lead- 
ing to the well-known Kubo formula: 



CTy(w) 



-y — 



1 



X (to| J,|n)(n| J,|to) (e"^^" - e"'^^'") . (C9) 



Comparison of Eq. ( C8 1 with Eq. ( C9 ) yields the general 



regularization procedure: 



E 



-I3E„ 



LOr, 



E 

(CIO) 

In a single-electron representation, the Gibbs factors 
are substituted for the Fermi occupation num- 



bers np{En). This procedure was used in Sec. IID to reg- 



ularize the EOM solutions of graphene in the presence of 
a magnetic field. 

In Sec. |II C[ no regularization was employed to derive 
the interband universal conductivity of graphene in zero 
field, Ke a^xi^j); see Eqs. (21) and (22 1. The reason 
is that the ^ pre-factor [coming from the electric field 
Ex{oj) = iloAq] is canceled by numerator in Eq. (CI) in 
this particular example, since in zero field, Jxi'^) ^ t^- 
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It is straightforward to show that applying the regular- stands, diverges. 



ization, Eq. (CIO I, to Eq. (21) yields exactly Eq. (22 1 



As for the imaginary part of the conductivity, the reg- 
ularization, Eq. ( |C10 I, is compulsory in order to obtain 
a consistent result; the imaginary part of Eq. (21 ), as it 



The regularization prescription, Eq. (CIO), is general 



and makes the link between the solutions of the EOM 
ijj{oj) [Eq. (C3l] and the exact regular optical conductiv- 
ity (t(cj) of electronic systems. 
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